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ABSTRACT 

The application of LORAN-C in the hyperbolic mode as a 
positioning system for hydrographic surveys was investi- 
gated. Observed LORAN-C time differences from a field test 
conducted in Monterey Bay, California were compared to 
calculated time differences determined from geographic posi- 
tions based on a microwave positioning system. Four methods 
were used to determine the calculated time differences. The 
first three methods were (1) applying only rhe seawater 
Secondary Factor, (2) computing the time difference based on 
a Semi-Smpirical TD Grid, and (3) applying ASF Correctors 
from the DMAHTC LORAN-C Correction Table. The final method 

applied multiple observed ASF Correctors at five minute 
latitude and longitude intervals. By applying mulriple 
observed ASF Correctors, which was the most accurate method, 
a 38.3 merer 1 drms with a lane offset of 3 to 12 merers 
using the 9940 X-Y LORAN-C combination was obtained. Based 
upon the results presented, it may be possible to use 
LORAN-C for hydrographic surveys at scales of 1:80,000. 
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I. 



I NTRO^CriON 



A. OSE OF LORAN-C 

In recent years there has been an increasing usage of a 
LORAN-C receiver and the LORAN-C network as the primary 
horizontal control for such scientific studies as deep ocean 
dumpsites, marine fisheries studies [Rulon, 1979], bathyme- 
tric surveys, and recently, a reconnaissance hydrographic 
survey. Examples of bathymetric surveys conducted by the 
National Oceanic and Atmospheric Administration (NOAA) , 
which have used Loran-C for positioning are: 

1) Su-100-1-79 Gulf of Alaska [NOAA H-9822, 1979], 

2) SO-100-2-79 Gulf of Alaska [NOAA H-9823, 1979], 

3) S-D902-WH-82 U.S, Nest Coast [NOAA Ship Surveyor, 

1982 ]. 

Recently an attempt was made to use LORAN-C as the 
sounding position control for a reconnaissance hydrographic 
survey S-K902-HH-82 [NOAA, 1982]. This was a special survey 
conducted by the NOAA Ship Whit ing in May 1982. Special 
surveys are field examinations of very limited extent or 
scope and frequently require unique survey or data collec- 
tion procedures [Umbach, 1976]. The purpose of this project 
was to verify the existence and extent of reported shoaling 
in three safety fairways in the Gulf of Mexico [NOAA, 1982]. 

The use of Loran-C as a positioning system for basic 
hydrographic surveys has been very limited due to the abso- 
lute accuracy of the long range system. A basic hydro- 
graphic survey is defined as a survey which is so complete 
that it need not be supplemented by other surveys. *'It must 
be adequate to supersede for charting purposes all prior 



9 



surveys" [Umbach, 1976]. Variables which affact rhe accu- 
racy of LORAN-C are signal propagation variations, weather, 
and sky waves. The effects of weather and sky waves on 
LORAN-C propagation are best described by Samaddar [1980] 
and the American Practical Navigator, [DMA, 1977] respec- 
tively. 

Signal propagation variations are due to the phase 
retardation of the signal as it passes over an all seawater 
path, over land paths, or partial seawater-land paths as 
compared to free space. Table I summarizes phase retarda- 
tion changes [Mortimer, 1 978]. Errors due to an all sea 

TABLE I 

Phase Retardation or Lag of Radio Waves 



1 ' 

Propagation Path 


Ren res entation 
Propa gation 
Vel ocit 7 
(kffl/sec) 


Difference in Phase 
Lag at 500 km 
Compared with 
Wave in Line 
Above (m) 


V acuum 


299792.5 


- 


Direct wave through 
earth's atmosphere 


299 691 


170 


Sround wave over 
sea water 


299 560 


220 


S round wave over 
rugged mountains 


298 899 


1,300 

1 



water path are known as the Secondary Factor (SF) and errors 
due to a land path or mixed path are known as the Additional 
Secondary Factor (ASF) [Speight, 1982]. 
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ASF Corrections in the LORAN-C system can be as large as 
plus or minus four microseconds, which is 600 meters on rhe 
baseline. In other areas with the same LOEAN-C coverage, 
these corrections may be much larger due to the expansion of 
the distance between adjacent hyperbolic lines of position. 
For example, at 32<> N and 8 0o W, using lattice patr 9960-X, 
a four microsecond (^sec) error will offset the 9960-X line 
of position approximately 2438 meters [Speight, 1982]. 

B. APPLICATION OF ASF CORRECTORS FOR NAVIGATION 

To compensate for the Loran-C positional errors caused 
by the ASF Correctors, the Secretary of Transportation 
tasked the Defense Mapping Agency Hydrographic/Topographic 
Center (DMAHTC) , the National Oceanic and Atmospheric 
Administration (NOAA) , and the United States Coast Guard 

(USCG) with the job of determining and applying the 

Additional Secondary Correctors for each Loran-C chain. 
This task was published in the Department of Transportation 
(DOT) National Plan for Navigation in the July 19, 1974 

Federal Register. These corrections should provide 95"5 

assurance that a vessel could fix its position to a pred- 
icted accuracy of 1/4 nautical mile (NM) within the U. S. 
Coastal Confluence Zone (CC Z) and the Great Lakes. The CCZ 
is defined as: 

"the seaward approaches to land, the inner boundary of 
which IS the narbor entrance and the cuter boundary of 
which is 50 nautical miles offshore or the edge of the 
Continental Shelf (100 fathom contour) whichever is 
greater. " 

The 1/4-NH accuracy requirement also affects the 

nautical chart. The National Ocean Survey (NOS) , which 
publishes charts for the CCZ, engaged in a program with the 
USCG and DMAHTC to provide the coastal navigator with charts 
overprinted with lattices which meet 1/4 -nm accuracy. The 
USCG, as operator of the LORAN-C radionavigation system. 
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conducts surveys to ensure that LORAN-C coverage exists 
within the CCZ and will be reponsible for the verification 
of 1/4-NM accuracy for all coastal LORAN-C service. In 
conjunction with NOS, it assists in surveys of coastal 
waters of the United States to allow producxion of LORAN-C 
charts based on observed field data to meet the standards 
set forth above (Speight, 1982]. 

DMAHTC, for LORAN-C civil need, prepares grid predic- 
tions from its data base. Based on analysis and verifica- 
tion of the predicted grid from a USCG and/or NOS survey, it 
produces revisions to the initial grid predictions [Speight, 
1982]. At present, DMAHTC has provided NOS with ASF 
Corrected LORAN-C Lattices which are overprinted on the NOS 
Charts. Each chart with ASF Correctors applied contains one 
of the following notes: 

"The LORAN-C lines of oosition overprinted cn this chart 
have been prepared for use with groudwave signals ana 
are presently compensated only for theoretical' prooaga- 
tion delays, wnich have not yet been verified by 

observed data. Mariners are cautioned not to rely 

entirely on the lattices in inshore waters. Skywave 

corrections are not provided". 

or 

"The LORAN-C lines of position overprinted on this chart 
have been prepared for use with groundwave signals and 
are compensated with propagation delays computed from 
observed data. Mariners are cautioned not to rely 

entirely on the lattices in inshore waters. Skywave 

corrections are not provided" [Speight, 1982]. 

Presently, all of the NOS Charts of 1:80,000 to 
1:120,000 scale covering the east coast. Gulf coast, and 
Great Lakes show LORAN-C lattices that have been compensated 
for Additional Secondary Factors. Most of the lattices on 
these charts have been constructed from DMAHTC data tapes 
that provide adjusted LORAN-C readings for each rate at 
every five minutes of latitude and longitude. A few 
lattices were constructed using a single ASF Correction for 
the entire chart area. Five minute data tapes were not 
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furnished by DMAHTC for constructing lines of position for 
LORAN-C rates on the West Coast Charts. On these charts a 
single average ASF Correction was used to adjust each 
lattice [NOAA, Marine Chart Division, 1982]. 

In addition to supplying corrected LORAN-C lattices for 
nautical charts, DMAHTC prepares, distributes, and periodi- 
ccdly updares unclassified ASF LORAN-C Correction Tables 
[Speight, 1982]. The ASF Correction Tables are for preci- 
sion navigation, utilizing digital computers to convert 
LORAN-C time differences to geographic coordinates [Speight, 
1982]. Presently, the ASF correctors found in the LOEAN-C 
Correction Tables were determined using theoretical propaga- 
tion delays. ASF correctors listed in the tables are going 
to be updated with observed data and reprinted the first 
quarter of 1983 [Wallace, 1 982]. 

C. APPLICATION OF ASF TO HYDROGRAPHIC POSITIONING 

Schnsbele [1979] investigated the possibility of using 
Loran-C as an electronic positioning system for hydrographic 
surveying. He concluded that in Monterey Bay, California a 
single Additional Secondary Factor (ASF) applied to offshore 
lines of position gave a root mean square error (drms) of 66 
meters for the West Coast 9940 Y-W pair and a predicted 42 
meter drms error for 9940 X-Y rates. 

The 42 meter predicted drms is larger than Nelson's 
[General Electric Co., 1979] findings in San Francisco Bay. 
He demonstrated, in a dynamic mode, that the precision of 
LORAN-C was 60.8 meters 2 drms (30.4 meter 1 drms) with a 
worst case of 71.2 meters 2 drms (35.6 meter 1 drms) . A 
mean difference or offset between the measured time differ- 
ence and the calculated time difference was 34 nanoseconds 
for the 9940-X rate and one nanosecond for the 9940-Y rate. 
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He also obtained a precision of 38.0 meters 2 drms (19.0 
meters 1 drms) in the static mode. Nelson also srates, that 
the above precision is only achievable if the user has a 
LORAN-C receiver which has the performance capabilities of 
those used in the experiment. The LORAN-C receiver must 
have "comparable signal averaging rime, extra notch filters, 
and attenuarion of the signal" [General Electric Co., 1979]. 

D. OBJECTIVES 

The National Ocean Survey requires that hyperbolic 
control systems used for hydrographic surveying exhibit a 1 
drms of less than 0-5 millimeter ar the scale of rhe survey 
[Umbach, 1976]- Although this requirement is generally for 
2 mHz phase comparision systems, it can be inferred that it 
also applies to other hyperblic systems such as LORAN-C. 
The scale routinely used 'for coastal surveys is between 
1:40,000 and 1:80,000 [Umbach, 1976] yielding an allowable 
error of 20 to 40 meters nor including systematic errors. 
Schnebele [1979] concluded that hyperbolic LORAN-C, after 
applying a single ASF Corrector, is unsuitable for basic 
hydrographic surveying. 

Whether or not applying multiple Additional Secondary 
Factors (ASF) to LORAN-C lines of position will reduce the 
drms sufficiently to meet the accuracy standards set by the 
National Ocean Survey Hydrographic Manual will be ascer- 
tained in this study. The term multiple ASF Correctors 
refers to the application of more than one corrector to 
LORAN-C lines of position over a given area. The variable 
ASF Correctors result from varying delays of the electromag- 
netic wave as it propagates over different land segments. 

Three methods of applying multiple ASF Correctors were 
tested. The first method was the application of a 
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Semi-Empirical Time Difference Grid Calibration Model devel- 
oped by The Analytic Science Corporation [1979]. The spon- 
soring agency was the Qnited States Coast Guard. The 

Semi- Empirical Model applies Secondary Factors and 

Additional Secondary Factors for each geodetic position 
based on the distance over land, the distance over water, 
and the total distance using mean sea water and land conduc- 
tivities. 

The second method which was investigated applies ASF 
Correctors found in the DMAHTC LORAN-C Correction Tables 
[DHAHTC, 1981] to LORAN-C lines of positions. These correc- 
tors were derived from the ground conductivities which have 
been determined in the field by a Coast Guard calibration 
team [D.S. Naval Oceanographic Office, 1982]. 

Finally, a third method was pursued. ASF Correctors, 
which were determined by field observations, were applied to 
the LORAN-C lines of position. These were determined by 
computing the difference between the observed LORAN-C rates 
and the expected time difference which was calculated using 
four lines of position from a very accurate microwave posi- 
tioning system. These ASF Correctors were determined at 
five minute latitude and longitude intervals. 
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II. NATUEJ OF THE PROBLEM 



A. THE PRINCIPLES OF LORAM-C 

To understand the problems associated with LORAN-C when 
used during hydrographic operations, one must first under- 
stand its principles of operation. LORAN-C is a low 
frequency, pulsed signal, hyperbolic, radio navigation 
system, employing time difference measurements of signals 
received by the navigator from at least three ground tran- 
smitting stations [Speight, 1982]. The stations are 
comprised of a master transmitting station, two or more 
secondary transmitting stations which are strategically 
spaced several hundred miles apart and, if necessary, a 
System Area Moniter (SAM) Station [U.S. Coast Guard, 1974], 

System Area Monitor (SAM) stations associated rfith each 
LORAN-C chain apply differential-type corrections to the 
rates in real-time. SAM stations continuously monitor the 
signals from all transmitters in the chain. If the oberved 
time difference deviates by - mors than 0.05 jusec from the 
expected value, then the appropriate secondary adjusts its 
emission delay time in order to remove the error [Schnebele, 
1 97 9 ] . 

The master and at least two secondary stations are 
located such that the signals from the transmitting stations 
can be received throughout the desired coverage area. The 
master station is designated by the letter "M" and the 
secondary stations or slave stations are aesignatisd W, X, 
Y, or Z [.O.S. Coast Guard, 1974 ]. 

All stations transmit on the common 
kHz. Interference between transmitters 
time separation [Poppe, 1982]. 



the use of 



16 



frequency of 100 
is avoided through 
After the master 



station transmits a pulse, each secondary station delays its 
own transmission for a fixed time, called the secondary 
coding delay. This coding delay is synchronized through the 
use of cesium frequency standards at each station. The high 
stability and accuracy of these standards permit each 
station to derive its own time of transmission without 
reference to another station [DMA, 1977]. Secondary coding 
delays are predetermined by system propagation times and 
equipment characteristics [ Laurila, 1976 ]. 

The pulse from the master transmitter is distinguished 
from those of the secondaries through phase coding of the 
pulses. Phase coding refers to the inversion of the nega- 
tive and positive peaks of the sine wave comprising the 100 
kHz carrier portion of the pulse. The purpose of the phase 
coding is twofold: 

"First, it permits automatic discrimination between the 
master and the various secondary stations, thereby 
permitting all stations to be identified by their rela- 
tive timing with respect to the master" 



"Second, the phase coding arovides protection against 
excessively long skywave delays which would cause the 
late arrival of the preceeding pulse to coincide with 
the leading edge or groundwave portion of a oulse being 
tracked" [Poppe, 1982]. 

The signals are received by a mobile receiver where the 
differences in time of arrival of the master signal and 
various secondary signals are measured and displayed on the 
indicator portion of the LORAN-C set. The accuracy of this 
time difference is increased by phase comparision "of the 
synchronized 100 kHz carrier within the master and secondary 
pulses" [Laurila, 1976]. This measured time difference (TD 
- in microseconds) represents a hyperbolic Line of Position 
(LOP) [U.S. Coast Guard, 19 74 ]. The intersection of .two or 
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Figure 2.1 Hyperbolic Fix (From Coast Guard LORAN-C User 
Handbook, 1974) 



more LORAN-C LOP's defines the 
(Figure 2.1). When plotted cn a 
the resultant hyperbolic lines 
position r Speight, 1982]. 



position of the observer 
chart, the intersection of 
defines a geographical 
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B. PHASE LAG 



In a vacuum, the velocity of radiated energy from an antenna 
for LORAN-C is 299792.458 Icm/sec. Since radiated energy 
cannot be shown pictoriall y, the phase of the transmitted 
radiations is used. The lines of constant phase of the 
transmitted radiation are shown in Figure 2.2 by the curved 
lines labeled aa ' , bb ' , and other similar designations. 
They define the wave front as it proceeds outward from the 
antenna in all directions. The distance between each line 
of constant phase is one wavelength (x) [Admiralty Manual, 
1965]. 



X (meters) = 29 9792.4 58 kra/sec e frequency in kHz 




Figure 2.2 Transmitted Radiation 



The velocity of the 
temperature, pressure. 



radiated energy 
humidity, and 



in air depends on 
the nature of the 
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surface over which the transmissions pass. The retardation 
of a transmitted wave is known as phase lag. When low 
frequencies are employed, such as LOHAN-C at 100. kHz, the 
effects of change of temperature, pressure, and humidity are 
swamped by the effects caused by changes in the nature of 
the surface over which the transmissions are traveling 
[Admiralty Manual, 1965]. The General Electric Company, 
TEMPO division, conducted a LORAN-C Signal Analysis 
Experiment under the direction of the U.S. Coast Guard. 
This experiment was conducted along the O.S. West Coast. 
The General Electric Company recorded an overall change of 
108 nsec and 116 nsec time of arrival from the master and 
X-secondary stations respectively after a storm [Samaddar, 
1980]. If ASF Correctors are as large as two microseconds 
on the West Coast [DMAHTC, 1981] then the ASF Corrections 
are 20 times larger in magnitude than weather effects for 
the 9940 West Coast LORAN-C chain. 

The change in transmission rates or phase lag are a 
result of the amount of energy transferred from ^he tran- 
smitted radiation. This energy transfer depends on the 
absorption qualities (inversely related to conductivity) of 
the surface over which they are passing and their wavelength 
(or frequency) . The lower the conductivity and the longer 
the wavelength (or lower the frequency) the greater the 
transfer of energy, and vice versa. Seawater has a rela- 
tively high conductivity (5.0 mhos/meter). Land has a much 
lower conductivity, which varies from marsh (fairly high) to 
dry sand and rock (very low) [Admiralty Manual, 1965]. Two 
excellant papers that discuss the electrical properties of 
soil are those of Smith-Rose [ 1934 ] and Pressey, Ashwell, 
and Fowler [1956]. Smith-Rose [1934] found that the conduc- 
tivities for soil ranged from 0.18 mhos/meter for a grey 
clay with salt to .00001 mhos/meter for granite. 
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Figure 2.3 Phase Lag 



Phase lag is illustrated in Figure 2.3. The lines of 
constant phase , aa\ bb', and cc' become distorted as they 
progress along the sea surface. The dotted lines represent 
the position of the lines of constant phase in the absence 
of rhe sea surface. The wavelength (X) in meters, measured 
at heights of several wavelengths above the sea, remains 
about the same as the direct wave through the earth’s atmo- 
sphere at 299691 hm/sec j frequency in kHz. Near sea level 
the absorption of energy retards rhe progress of rhe wave- 
front, and makes the sea level wavelength (A') less than A. 
As the lines of constant phase progress away from, the 
antenna the phase lag increases wi-h distance. This is 
known as the Secondary Factor (SF) [Admiralty Manual, 1965]. 

The most intriguing fearure of phase lag occurs an the 
coastal interface where there is an extreme change in the 
con duct ivites between land and sea. Visualizing the wavs 
front in three dimensions, the lower part of the wave, 
slowed by the drag of the ground, lags further and further 
behind the upper par~ as the wave crosses the land. At the 
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coastline it suddenly encounters the much lower impedance of 
the sea, and in a very short distance the bottom of the wave 

tries to catch up with the top, as though the whole wave 

front were an elastic balloon. This is known as "phase 
recovery" [Eaton, 1979]. Phase recovery was verified during 
tests on Decca transmissions across the south coast of 
England by Pressey, Ashwell, and Fowler [1956]. 

The determination of the Secondary Factor for seawater 
is fairly direct since the conductivity of seawater (5.0 
mhos/meter) is fairly constant. But for land the conduc- 
tivity can vary depending on the type of soil and its water 

content [Smith-Rose, 1934]. 

Phase lag for radiated energy over land can be deter- 
mined two ways: 

1) Assign an average land conductivity to the ASF 
Model. For example, the average conductivity for 
the soil on the west coast is 0.003 mhos/meter. 

The average land conductivity will determine the 
average phase retardation of the path [The 

Analytic Science Corporation ,1979]. 



2) Determine every conductivity for each portion 



of a line segment from 
receiver. The total 
segments constituting a 
determine the total phase 
[Speight, 1982 ]. 



the transmitter to the 
of these conductivity 
land-water profile will 
retardation of the path 



C. TC MODEL 

Positional fix accuracy using LORAN-C is primarily 
dependent on a chart makers ability to accurately compute 
the expected difference in time-of-acri val (T0A> of received 
groundwave signals from the transmitting stations. Time 
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differences (TD) , are the differences between the TOAs of 

the secondary and master transmitters. 

TD = TOA - TOA (2.1) 

i i m 

i = Secondary Station 
m = Master Station 



TOA computations are dependent upon an accurate knowledge of 
the signal phase delay. 

The phase delay of a groundwave signal is generally 
expressed as: 



$ 



T + SF 



nR 

C 



+ SF 



( 2 . 2 ) 



where n is the surface refractive index, C is the speed of 
light in a vacuum, R is the range between the receiver posi- 
tion and the transmit ting station. The primary phase delay, 
T, is the computed travel time of the LORAN-C pulse over a 
distance equal to the tran smit ter-to-receiver greax. circle 
path length, taking into account the velocity of electromag- 



netic waves and the index of : 
The secondary factor (SF) is 

phase delay and accounxs for 
term in (2.2) is the primary 
usually an order of magnitude 
Corporation, 1979 ]. 

Thus, time-of-arrivals can 

TOA = T + 

i i 

TOA = T + 

m m 

where CD is the xrue emission 
LORAN-C chain [The Analytic Sci 
coding delay is equal to a 



efraction of the atmosphere. 

a correction to the primary 
he phase lag. The dominant 
phase delay (T) . The SF is 
smaller [The Analytic Science 

be expressed as: 

SF. + CD. (2.3) 

SF (2.4) 

m 

delay or coding delay for the 
ence Corporation, 1979]. The 
me delay plus a compuxed one 
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way baseline time (Be) which includes the secondary phase 
correction for an all seawater path. The oneway baseline 
time (Be) is equal to the distance between the masrer and 
secondary transmitters in meters divided by the propagation 
velocity of LORAN-C through the earth's atmosphere (299.691 
meters per microsecond [Navigation Department DMA, 1982]). 
See Table II for Ceding Delay values for the 9940 chain 

TABLE II 

Coding Delay 9940 LORAN-C Chain 



Pair 9940-W: CD + 3c = 11000 ♦ 2796.90 = 13796.90 usee j 

9940-X: CD + Be = 27000 + 1094.49 = 28094.49 jusec | 

9940-Y: CD + Be = 40000 + 1967.27 = 41967,27 usee | 

I j 



[Riordan, 1979]. Combining equations 2.1, 2.3, and 2.4, the 
true TD is given by equation 2.5 [The Analytic Science 
Corporation, 1979 ]. 

TD = (T - T ) + (SF - SF ) + CD (2.5) 

i i m i m i 

D. SEMI -EMPIRICAL TD GRID CALIBRATION MODEL 

The Semi-Empirical TD Model was developed by The 
Analytic Science Corporation [1979] in Reading, 
Massachusetts for the Sest Coast. 9940 LORAN-C chain. 
Similar "time difference (TD) grid calibration techniques 
have been successfully employed to develop an accurate 
(approximately 100 nsec drms) calibrated grid for St. Marys 
River LORAN-C chain", [The Analytic Science Corporation, 
1979 ]. 
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1 • Te chnical A ppro a ch 



The Semi-Empirical Model is based on Millington's 
empirical approach for computing the secondary factor over a 
mixed (multiple-homogeneous segment) path which combines 
land and sea phase delays. The generalized semi-empirical 
polynomial functional form for the SF of the LOHAN-C station 
is given by: 

K 

2 L 

SF = SF(T ,6) = Z AT * Z (C SinlS + D ccsl 6 ) (2.6) 

j j D k j 1=1 jl j jl j 

where 

j = secondary (WfX, or Y) or master (M) station, 
nH . 

T = — - = jth station-to user primary phase delay, 

j c 

R = jth stat ion-to-user great-circle path length, 

3 

B = user path bearing angle at the jth station, 

j 

K , K and L are positive integers, 

12 

C and D are the station-dependent coefficients 

jl 

of harmonic terms in the model, 

A is the range-dependent coefficient of the model 
k y t' 

which may in general be station-dependent. 

Data from 27 coastal sites distributed along the 
West Coast and 122 land-sea sites distributed in the 
Southern California CCZ (between Point Arguello and San 
Diego - see Figure 2.4) were used in a Kalman estimation 
algorithm to compute the uncertain coefficients of the land 
and sea models of the ID grid calibration algorithm. (An 
explanation of Kalman filtering for the layman is presented 
by Roger M. du Plessis [1967].) The calibrated algorithm 
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was used to compute TDs at each data site and the TD resi- 
duals (difference between measured and calibrated TDs) were 
examined. Adjustments were then made to the TD model struc- 
ture in an attempt to further reduce the residuals. This 
process of adjusting the model structure is repeated until 
rhe residuals agree with the expected theoretical covariance 
associated with the TD model. The model which exhibited rhe 
"best" performance was selected as the West Coast TD grid 
model. 

Gen eraliz ed R an g e/ B earing Model 

The Generalized Range/Bearing (GRB) Model was 
selected as the "best" semi-empirical calibration model for 
the West Coast chain. The semi-empirical function is: 

L 

S? = AO + AIT + z (C sinlS + D cosl 8 ) usee (2.7) 

j j 1=1 jl j jl j 

where AO, A1, C and D are the model coefficients, is the 
path bearing angle measured positive clockwise from north at 
the jth ( W,X,Y or M) station and T is the path range to the 

jth station. The GRB model is relatively complex and is 

expected to exhibit superior performance. The extensive 
model is based on knowing the distance overland (TL) , the 
distance over water (TS) , the total distance (T) , and the 
path bearing angle 8 (Figifre 2.5). 

It was noted that the calibrated model was expected 
to be accurate and applicable only over the extent of ranges 
and bearing angles embodied in the calibration data. Hence, 
outside the region covered by the calibration data the model 

may not be as accurate as within the data coverage region. 

Using the GRB model, a drms value of approximately 0.8 ^sec 
was expected in areas where land data alone was used to 
calibrate the model. Inclusion of sea calibration data 
produced a drms value of 0.35 to 0.50 psec. 
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NORTH 

I 




Figure 2.5 Nixed Path TD Geonserry 



The time difference (TD) is expressed by rhe 
following equation; 

TD = (T - T ) + (SF. - SF ) + CD. + b. /isec (2.8) 

i i m i m r i 
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whe re 



n R 

i 

T |isec, 

i c 

n R 

m 

T ;isec, 

m c 

R, = ith secondary srat ion-to-user great-circle path length. 



c = rhe speed of an electromagnetic wave in a vacuum 
= 2.99792458 X iQs m/sec, 

n = surface refractive index 
= 1.000333, 

R = master station-to- user great-circle path length, 
m 

CD.= ceding delay found in Table II, 



b , = TD bias associated with the ith secondary 
station (jasec) (Table III) . 



TABLE III 



TD Bias (b) - /isec 



TDW -0.854 1 

TDX -1.173 I 

TDY -0.353 1 

c. - - I 



SF.= 0.5 (-S1 + S2 + S3 - 34 + S5 + S6) (2.9) 



The term Si is -he SF of a land path of length T 

s j 

ijisec) from the jth station: 



Si = 0.795/T 

sj 


+ 0.439 


+ (0.00245) T 

s 1 


/isec. 


(2.10) 


The terms 


S2, S3 


, and S4 combine 


to make 


up the 


secondary 


factor 


for the seawater 


path le 


ng- hs. 
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S3 is rhe SF 



S2 is the SF for the toral path, 
using the seawater coefficients for the portion 
with seawater, and S4 is the SF using the seawater 
coefficients for the land path distance. The 
seawater coefficients are found in Table IV. 



TABLE IV 

Seawater Coefficients 



a1 = 128.8 
a2 = 0.187 

a3 = 0.000652 

b1 = 3.188 

b2 = -0.594 

b3 = 0.000329 

I I 



SF (T) = 

c; 


bl/T + b2 + (b3) 


T 


/IS 




if 10 < T < 540 




/isec 


SF (T) = 

o 


al/T + a2 + (a3) 


T 


;is 




if T > 540 ^sec. 







(2.11a) 
(2.1 1b) 



Term S5 and S6 are the SFs of land paths of length 
T and T . SF5 is the Secondary Factor for the 

j 1 j 

total lengrh using the land coefficients whereas 
SF6 is the Secondary Factor for the distance over 
land using the land coefficients. The land 
coefficients are found in Table V. 



55 = SF (T , 6 ) , 

L j j 

56 = SF^ (T^ ..6.) , 

L Ld :) 

whe re : 

SF (T , 3 ) = AO + (A1) 
L 3 j 



2 

T + E (C sinl 3 + D cosl 8 ) (2.12) 

j 1=1 jl j jl j 



30 



TABLE V 



Land Coefficients 



AO 


= 


1. 428 


D 

x2 


= 


0. 942 




A1 




0. 00158 


C 

y2 




o 

• 

o 




C 

w1 


= 


O 

• 

o 


= 


0.588 




c 

w2 




-0. 71 1 


D 

yi 


— 


o 

• 

o 




D 

Wl 




0. 323 


D 

y2 


— 


o 

• 

o 




D 

w2 




o 

• 

o 


c 

m 1 




1.010 




C 

Xl 




o 

• 

o 


C 

m2 


— 


-0.196 




c 

x2 




0. 0 


D 

m 1 


— 


-0.893 




D 

Xl 

1 




o 

• 

o 


D 

m2 




-0.355 


J 



E. DMAHTC MODEL 

1 . S§ a S F Mo del 

The equations for the Sea SF Model is: 

SF = (B1/T) + B2 + (B3 T) usee, if ID < T < 537 ;isec, (2.13a) 
SF = (A1/T) + A2 + (A3 T) u sec, if T > 537 ;isec, (2.13b) 

where T is rhe primary phase delay (or range) in microse- 
conds (usee); AJe and Bic (k = 1 , 2, and 3) are -he sea model 
coefficients used by DMAHTC in program TDGRID [Funakoshi, 
1982]. The coefficients are found in Table VI, 

2 . La nd SF Model 

The solution used to resol7e DMAHTC Tables ASF 
Corrections is called Millington's Method [DMAHTC, 1981]. 
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TABLE VI 



Sea SF Model Coeffioients 



A 1 




129.04323 


A2 


= 


-0.40758 


A3 


= 


0.00064576813 


B1 


— 


2. 741282 


B2 


= 


- 0. 011402 


B3 




0.00032774815 



I I 



This method is based on the premise that the phase distor- 
tion due to a composite land-sea path is the arithmetic 
average of the phase distortion found in the forward and 
reverse paths of the propagated signal [ DMAHTC, 1981]. For 
example, in Figure 2.6 two azimuths have been drawn on the 
map and are labeled as 210° and 235o. Also placed on the 
map are the proper ground conductivities which have been 
determined in the field by the Coast Guard calibration team. 
A great circle drawn on the appropriate chart or charts from 
the LORAN-C Station coordinates to the area under considera- 
tion spans various lengths of land and seawater. Each 

length or segment will have a specific conductivity and 
distance. The total of these conductivity segments, consti- 
tuting a land-water profile, will determine the total phase 
retardation along that path [U.S. Naval Oceanographic 
Office, 1982]. All azimuths and distances are computed 
based on the World Geodetic System (WGS) datum [DKAHTC, 
1981]. The values of phase retardation for a given ground 
conductivity are tabulated in the National Bureau of 
Standards (N3S) Circular 573 [Speight, 1982]. 

The formula used to derive the ASF Correction for 
the time difference for a master-slave transmitting station 
pair is: 

ASF Correction = (-Slave Error) - (-Master Error) 
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1 



TRANSMITTER 




Figure 2.6 Composite Laad-Sea Path 



The 

forward and 
subtracted 



mean values derived for one station from 
reverse solution of Millington's Method 
from the Sea SF Model. The differences 



the 

are 

are 
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presented in an azimuth array. This array is a series of 
geodetic azimuths radiating from the transmitter with 
corrections computed at incrementea distances along each 
azimuth. After the values are computed in the azimuthal 
array a computer program rearranges them into a matrix form. 
The matrix form is the arrangement of corrections into rows 
and columns covering a specified geographical area at a 
constant spacing. The purpose of the matrix is to enable 
the corrections from two LORAN-C transmitters to be added 
algebraically, combined into a single matrix, and arranged 
in the desired tabular form. This tabular form is the body 
of the table [DMAHTC, 1981]. 

3 • Ta b le D escripti on 

Each table contains a complete chain. Figures 2.5, 
2.6, and 2.7 depict LOHAN-C ASF Correctors for chains 
9940-H, X, and Y for Monterey Bay, California. A table 
section is prepared for each station pair (master station 
and one slave station) in a LORAN-C chain. As a rule the 
limits of the table coverage are determined by the range of 
the ground wave transmissions for the LORAN-C chain. Each 
page of corrections in the table covers an area three 
degrees in latitude by one degree of Icngirude, with correc- 
tions printed in increments of five minutes of arc. Rate 
designation and page numbers are printed at the top of each 
correction page. Those pages where latitude and longirude 
limits contain both land and sea are included bur correc- 
tions apply only fcr the area covered by the CJ.S. Coastal 
Confluence Zone (CCZ) . Large land bodies and areas ouside 
the CCZ are represented by blank spaces on the page. ASF 
correction values can be either positive or negative (posi- 
tive values are shown without sign) . Areas requiring no 
correction show a zero value which in some cases in preceded 



34 





LONGITUDE WEST j 

0 o 

122 121 i 

0' 55 50 45 40 35 30 25 20 15 10 5 0 ' 


9 < 




39 0 




55 




50 




«5 




40 




35 




30 




25 




zo 


1 


15 




10 




L 0 5, 




A 38 0 




T 




1 55 




T 50 




U 45 




0 40 




E 35 




30 




25 




20 




N 15 




0 10 


! 


R , 5 


1 

1 


T 37 0' 


1 

1 


H 


1 


55 


- 1 . 6 - 1 . 6 1 


50 


-1 .6 -1 .5 -1 .4 1 


I ^5 


-1 .5 -1 .4 -1 .6 \ 


1 40 


-1 .4 -1 . 3 -1 .5 i 


1 35 


-1.3 1 


1 30 


-1.4 i 

I 


1 25 


-1 .4 -1 .5 i 


1 20 


-1 .2 -1 . 1 ^ 


15 


-1.4 -1.3 -1.6 


10 


-1.3 -1.3 -1.7 -1.7 


0 5 


-1.3 -1.4 -1.7 -1.7 -1.8 


36 0' 


-1.2 -1.6 -1.6 -1.7 -1.7 -1.7 -1.8 






Figure 2.7 9940-W A3F Correcr-ion Table 



by a negative sign indicating that the zero results from the 
rounding off of a value slightly less than zero (indicates 
the trend of the correction). 



35 




I 








ship' s 
arc in 



The table can be entered directly by 
position determined to the nearest five 
latitude and longirude either by dead 



using the 
minutes of 
eckoning or 
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Figure 2.9 9940-Y AS? Correction Table 



some other 
correction , 
utilized . 



means. To find the page with the appropriate 
the Page Indexes of the table should be 
These indexes show the limios and page number of 
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all pages in the table. To locate the number of the page on 
which the desired correction is to be found the Page Index 
is entered with the ship' s position. In some cases the 
ship's position will fall on the page limit in either lati- 
tude or longitude or both. These positions are repeated on 
both pages and either page may be used. 

The ASF Correction is added algebraically to the 
time difference for the LORAN-C pair. Interpolation of this 
data will not necessarily improve the accuracy due to the 
method used to determine ASF Correctors [DMAHTC, 1981]. 
Since the correctors are computed in the azimuthal array and 
are based on the conductivity and distance over which the 
LORAN-C electromagnetic wave travels, the ASF Corrector 
between the published ASF Correctors in the tables may not 
be the linear interpolated values. For example, the ASF 
Corrector for a distance of 500 m with an azimuth of 180° is 
equal to 1.5 jisec. The ASF corrector for a distance of 500 
m with an azimuth of 181<3 is equal to 1.6 psec. The inter- 
polated value between 1-5 and 1.6 is 1.55. The true ASF 
Corrector for the 180° 30' azimuth is 1.4 since the land 
distance for the same azimuth is less than the land distance 
for the 1800 and 181° azimuth. The LORAN-C signal passed 
over a harbor [Dansford, 19 82 ]- 

F. ATTEMPTED DETERMINATION OF ASF CORRECTORS BY HYDRO FIELD 
PARTIES 

One of the major problems encountered by hydrographic 
survey operating units when using LORAN-C for position 
control is the determination of the AS? Correctors for the 
survey area. The four surveys mentioned in the Introduction 
all made attempts to determine the correctors by comparing 
the LORAN-C rates to a second source. 
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Bathymetric Surveys H-9822 [NOAA H-9822, 1979] and 

H-9823 [NOAA H-9823, 1979] Gulf of Alaska, compared the 

rates from an Intern av LC-204 LORAN receiver to computed 
rates from a position obtained from a JMR-1 Satellite 
Navigation Receiver whan available. Shore ties using radar 
ranges, visual bearings, and sextant angles in comparison to 
LORAN-C rates were also made prior to and after each survey. 
The calibrations of LORAN-C rates were based on the satel- 
lite positions only since the positional computation of 
LORAN-C and JMR Doppler Satellite were made on the WGS 1972 
datum whereas the land ties were based on the NAD 1927 
dat urn . 

Bathymetric survey SU- 40-7-82 which extended along the 
Washington, Oregon, and California Coasts used LORAN-C as 
navigational control. LORAN-C time differences were 

compared with SATNAV positions. The report did nor indicate 
whether any correctors were applied [NOAA Ship Sur vey or, 
1 98 2 ] . 

Finally, Hydrographic Survey S-K902-Wh-82, 

Reconnaissance Survey of Safety Fairways, Gulf of Mexico 
used LORAN-C as a posirioning control system. The posi- 
tioning unit was an LC-204 receiver. LORAN-C rates were 
input via the HIDROPLOT Controller, a special purpose 
input-output interface which is the nucleus of the computer 
sysrem hardware [Ombach, 1976]. Positions were computed and 
plorted by Program P.K121, LORAN-C Real-Time HYDROPLOT 
[ Backus , 1 980 ]. 

ASF Correctors for LORAN-C were ichieved by visual cali- 
bration using three poinn sextant fixes using charted oil 
rigs as control in the vicinity of nhe survey area. A 
three-point sextant fix is a convenient and accurate method 
for determining the position of a hydrographic survey 
vessel. Sextants are used to measure two angles between 
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three objects of known geographic positon. The center 
object is common to both angles. The position of the obser- 
vers taking the angles is fixed by the intersection of three 
circular lines of position [Umbach, 1976]. 

These sextant angles ware recorded and later tranferred 
onto their respective charts using a plastic three-arm 
protractor. A plastic three-arm protractor is transparent 
and made up of one fixed arm and two movable arms which 
contain an etched line that is radial with the center of the 
protractor [Ombach, 1976]. Sextant angles observed in pairs 
for a resection fix with a common center mark may be plotted 
directly by this instrument. When the three arms are placed 
at the angles observed and fitted so as to pass through the 
plotted positions of the observed stations on the field 
sheet, the hole at the center of the three-arm protractor is 
the fixed positon of the vessel [Ingham, 1975]. 

Partial correctors for each area surveyed were defined 
by comparing the observed rates and the determined rates 
plotted on the nautical chart. The partial correctors were 
applied via the HYDRO PLOT Controller. However, even after 
applying these correctors, the plotted position still disa- 
greed with the ship's determined pcstion with respect to the 
oil rigs. Ship's personnel attributed the discrepancies to 
one or more of the following: 

1) Accuracy of the charted rigs, 

2) Weather effect on LORAN-C, 

3) Time of day, 

4) Propagation of signal over land path, 

5) Three-arm protractor accuracy, and 

6) Error in the conversion by the software 
of the LORAN-C rates to latitude and 
longitude [NOAA, 1982]. 
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There is an apparent need for a LORAN-C calibration 
routine aboard NOAA ships which provides the ASF Correctors 
for program RK121, LORAN-C Real-Time HYDROPLOT. The routine 
should use the same geodetic distance computation found in 
RK121 and use the same datum as that of the nautical chart 
of the survey area. 

The above mentioned discrepancies illustrate the defi- 
ciencies in applying a single ASF Corrector zo LORAN-C data. 
The accuracies for hydrography cannot be met using single 
correctors because the errors are non-linear and systematic. 
They cannot be distributed like residuals in a traverse. 
Schnebele [1979] has already proven that single ASF area 
correctors to LORAN-C positions do not meet the accuracy 
standards of the NOS Hydrographic Manual. 

Based on visual inspection of the DMAHTC LORAN-C 
Correction Tables, ASF Correctors should be updated every 
five minutes of Latitude or Longitude change. In Monterey 
Bay, California, there is approximately 0.1 to 0.2 ^sec 
difference for every five minutes of change, a potential 
error of 55 to 110 meters. 
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III. EXPERIMENTAL PROC EDU RE 
A. FIELD PROCEDURES 

In order to compare the use of differential LOHAN-C with 
ASF multiple correctors, typical survey operations were 
planned for the southern portion of Monterey Bay, 
California. This survey was conducted in conjunction with a 
comparative evaluation of multiple lines of positon for 
selective positioning methods [Anderson, 1982]. Four micro- 
wave ranging systems were set on known geographic positions 

TABLE VII 

Geographic Names and Positions 



Microwave System 


Ge oarapi 


lie Position 




; used. 


Stations 


(NAD 1< 


327) 


Dates 


Seaside 4 <1 964) 


360 


36 • 


23.44596” 


June 


3-5,1982 


121 0 


51 ' 


38.83281” 






Use Mon Ecc. 


360 


36 » 


04.73031” 


June 


3-5,1982 




121 0 


52’ 


35.98040” 






Geoceiver See. 


360 


36 ’ 


32.49281” 


June 


3-5, 1932 




121 0 


53’ 


25.21162” 






Mussel Ecc. 


36 0 


36 ’ 


18.25484” 


June 


3-5, 1982 




1210 


54’ 


11.49661” 






Park (1931) 


36 0 


53’ 


13.80600” 


June 


6-7,1982 


121 0 


49’ 


46.74300” 






Mulligan RM1 


36 0 


44 ’ 


56.49531” 


June 


6-7,1982 


121 0 


47’ 


52.31090” 






Range 7 (197 2) 


36 0 


39’ 


02. 47787” 


June 


6-7, 1982 


121 0 


49’ 


08.58202” 






Mussel (1 932) 


36 0 


37’ 


18.15100” 


June 


6-7,1982 


1210 


54’ 


11.49661” 




i 
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listed in Table VII. A series of tracklines were run in two 
separate areas as shown on Figure 3.1. To ensure that the 
microwave postioning system was working properly, the equip- 
ment was calibrated over known baselines of 1497.47 meters 
and 7877.31 meters at the beginning and end of the project. 
Trackline observations were only made during the daytime in 
fair weather conditions so as to eliminate sky waves and 
weather changes that influence L9RAN-C signal propagation 
characteristics [Samaddar, 1982]. The vessel used was the 
126 foot R/V Acania which is operated by the Naval 
Postgraduate School. 

The positioning equipment consisted of a Micrologic 
ML- 1000 LORAN-C receiver (0.01 usee resolution) and a 
Trisponder Microwave System provided by Racal-DECCA Survey, 
Inc. The Trisponder Microwave System consisted of four DNT1 
Model 217C transponders, four DNT1 Model 21017 HP sector 
antennas with 870 by 50 beam widths, one DNT1 DDMU (Digital 
Distance Measuring Unit), two Omni DVTI Model 21019 HP 
antennas, a Houston Instruments Model DP3-M2D/RC3 plotter, 
and a Texas Instruments 7 43 terminal (Table ’/III). The 
manufacturer's published accuracy for the positioning equip- 
ment is ± 1 m for a single range [Racal-DECCA Survey, 1981]. 
Anderson [1982] discusses the accuracy of four lines of 
position. The four Decca Trisponder distances were recorded 
via a Texas Instruments 743 data terminal while the LORAN-C 
rates were manually logged . The data was acquired at one 
minute intervals while the ship maintained constant course 
and speed. The recorded LORAN-C rates were 9940-W, X, and Y 
of the West Coast chain. 

To test the potential for calibrating the LORAN-C System 
using the Semi-Empirical Model, the correction tables, and 
multiple observed field correctors, the positions derived 
from the microwave system measurements were used to compute 
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I 




TABLE VIII 



Microwave Positioning Equipment 



E quip me nt 

DNT1 Model 217C Transponders 



(870 by 50 Beam Widths) 

DNT1 DDMV Model 540 
Omni, DNT1 Antenna Model 2U 

Houston Instruments Plotter Model DP3-M2D/RC3 
Texas Instruments 743 Terminal 





S/N 


Code 72R] 


3323 


Code 74R; 


3320 


Code 76R 


3321 


Code 78R] 


3322 


nna 


185 




186 




187 




191 




426 


HP 


194 




200 



10722-10 

34418 



expected LORAN-C time differences ar each point. The 
difference or offset between these expected time differences 
and the observed values were computed for the nhree merhods. 
The mean offset, standard error, and drms values were also 
computed and compared. 

3. MICROWAVE SYSTEM POSITIONING 

The geographic position of the ship based on four lines 
of position was determined using a computer program called 
GP3YLQ (Geographic Posirion by Lease Squares) written by the 
author (See Appendix A) . GPBYLQ contains subroutine LSQR 
(Least Squares) , which is a least squares adjustment written 
by Paul R. Wolf, Ph.D. [1974] and revised by LCDR D. Leath 
[1981]. Geographic postions were converred to X,Y (meters) 
which in turn were converted to geographic position via 
subroutines GPTOXY and XYTOGP, respectively [Wallace, 1974]. 
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Subroutine GPTOXY and XYTOGP are based on the Modified 
Transverse Mercator Grid which was centered in the survey 
area. 

The Modified Transverse Mercator (MTM) projection is 
used by the National Ocean Survey and is similar to the 
projection used in the Universal Transverse Mercaror (UTM) 
system. The main difference is that in the MTM a Central 
Meridian is picked that is near the survey area instead of 
being fixed at a particular meridian [Wallace, 1971]. 
Central Meridian (CMER) , False Easting (FEST) , and 
Controlling Latitude (CLAT) are the three parameters which 
define the MTM projection. CMER is the mean longitude 
computed using the maximum and minimum longitudes of the 
survey limits, FEST is the X-Coordinate that is assigned to 
the Central Meridian, and CLAT is the distance in meters 
from the equator to seme reference latitude [Wallace, 1971]. 
The Central Meridian, False Easting, and Controlling 
Latitude used for Monterey Bay, California referenced to SOS 
Chart 18685 are; 

CMER = 121 0 56» 00.0", 

FEST = 20000.0, 

CLAT = 4050000.0, 

To be consistent with the National Ocean Survey charts 
of the area, all computations were done relative to North 
American Datum (NAD) 1927 geographic positions. All 
programs were executed on an IBM 3033 computer located at 
W.R. Church Computer Center, Naval Postgraduate School, 
Monterey, California. 

C. LORAN-C COMPUTATIONS 

The differences or offset between the observed and 
computed LORAN-C rates using the Semi-Empirical TD Model, 
ASF LORAN-C Correction Tables, or the Multiple Observed 
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TABLE IX 



LORAN-C 9940 Chain Data 



Geographic Position 

Staton (NAD 1927) 



Master - 


Fallon, Nevada 


390 

1180 


33' 

49' 


07.03»N 

52.23"W 


Slave - 


George, Washington 
9940- W 


O O 


03* 

44* 


48.82"N 

34.78»W 


Slave - 


Middletown, 
California 99 40-X 


380 

1220 


46 ' 
29» 


57.49»N 

40.04»W 


Slave - 


Searchlight, Nevada 
9940-Y 


350 

1140 


19 * 
48* 


18.32"N 
13.95* W 



Correctors were compared to the offsets between the observed 
TD rates and the computed rates for which only the seawater 
Secondary Factors (SF) were applied. The comparison of the 
offsets between the four methods illustrates the improvement 
in positional accuracy after applying ASF Correctors. 

1 . Sea^^r Sec on dary Fact ors (SF) 

Time differences using only the seawater Secondary 
Factors for each of the geographic positions were computed 
using program LORAN written by the author (Appendix C) . 
Seawater Secondary Factors (SF) were computed using formula 
2.13 and the coefficents found in Table VI. All TOA 
distances in meters were determined using subroutine INVERT. 
INVERT is a geodetic inverse routine using T. Vincenty's 
modified Rainsford's method with Helmert's elliptical terms, 
programed by LCDR L. Pfeifer, NOAA [1975]. Subroutine 
INVERT is accurate to 0.0001 m halfway around the world 
[Pfeifer, 1982]. All distances were converted to microse- 

conds using 299.792458 m/fisec. Time differences (TD) were 
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computed from equation 2.5. North American Datum 1927 
geographic positions were used for all computations. See 
Table IX for the positions of LORAN-C 9940 transmitters 
[Riordan, 1979]. 



To ensure that subroutine INVER1 was functioning 
properly the distances between the master and secondary 
stations were compared to the NOS published baseline 
distance [ Riordan, 1979]. The published distances and rhe 
results from routine INVER1 are lisred below. 



9940-W 

9940-X 

9940-1 



Publi shed B asel ine 
Dista nce (m) 
837,777. 0929 
327,886. 3720 
589,298. 5712 



Commuted Base line 
Distance - INVERT, (m) 
837,777. 115 
327,886. 316 
589,298. 589 



The difference between the published and computed baselines 
ranged from 0.02 to 0.06 m. 

Differences or offsets (x ^ ) were obtained by 
subtracting the observed LORAN-C rates from the compured 
values from the various methods. The mean difference or 
offset (X) and standard deviation (s) in microseconds for 
each rate were determined using equations 3-1 and 3.2 
[Wonnacott, 1935]; 



n 

1 

Mean offset (;asec) : x = - S x (3.1) 

n i 

i=1 

where; x = original observation in usee, 
i 

n = number of observations; 



Standard 

error: 

(>isec) 



n 

1 

s = E (X. - X) 2. (3.2) 

n- 1 1 

i = 1 
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The mean offset in microseconds can be converted to meters 



using equation 3.3 [Heinzen, 1977]; 



Mean offset 
(meters) : 



1 



X w 
sin a 



(3.3) 



where w = distance corresponding to one microsecond 

on the baseline 

= 149.896229 meters [Bigelow, 1963], 

a = one-half the angle between the 

radius vecrors from the position to 
the master and secondary stations. 



Equation 3.4 was used to compute drras values in meters 
[Heinzen, 1977]; 



drms= 



sin ( + ) 

1 2 



( s w) 2 (S W) 2 2 pcos ( a + a ) s w s w 
1 2 12 12 
+ + ( 3 . 4 ) 

2 2 (sin a ) (sin a ) 

sin a sin a 1 2 

\1 1 2 



where: p = correlation coefficient = 0.33, 



s and 



a are as above with the subscript 
denoting the appropriate 
secondary station. 



The correlation coefficient (p) is a result of the secondary 
station having a common line of position with the master 
station. Although often ignored, various authors assign 
values ranging from 0.33 to 0.40. Bigelow [1963] chooses 
P = 0.33. 



2 . Se mi-Empirica 1 TD G rid 

To determine if the Semi-Empirical TD Grid would 
reduce either the drms value or the offset between the 
observed and expected TD rates, program LOPLC (Line of 
Position - LORAN-C) was written by the author (Appendix B) 1 
Program LOPLC computes the distance over land, distance over 
water, the total distance, and the azimuth from north for 
the transmitting station using Subroutine INVER1. The 
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land-sea distances were computed by selecting a point which 
was located along the coast of Monterey Bay. Subroutine 
SHORPT (Shore Point) interpolated a geodetic position from 
23 geodetic points which outlined Monterey Bay, California. 
The land/sea Secondary Factor was determined in Subroutine 
SECFAC which is found in program LOPLC using equations 2.8 
through 2.12. All compu tarions ware based on NAD 1927 
geodetic datum. 

Tests were made on program LOPLC using data found in 
the Semi-Empirical TD Grid article [The Analytic Science 
Corporation, 1979] using the HGS 1972 datum. Station TASC 
55 located at latitude 34° 34' 18.3” N and longitude 120<> 

39' 40.3” W, was selected from The Analytic Science 

Corporation article. It was one of the stations used to 
calibrate the coefficients for the Semi-Empirical TD equa- 
tions discussed earlier. The only distances listed for TASC 
55 were the individual distances over land and over water. 
The total distance between the transmitters and TASC 55 was 
computed by adding the land and sea distances. The 
following station-tc- site path segment lengths are listed 

INVERT 
Tot al 
Distance 

(isa) 

576 . 083 
495.178 
541 .308 

Using the land distances, sea distance, and total distances 
to TASC 55, and the computed azimth from subroutine INVERT, 
program LOPLC produced offsets of 0.86 psec and 1.13 ^sec 
for 9940-X and 9940-Y respectively when compared to the 
expected time differences at TASC 55. 



for TASC 55: 



Station 
Master 
X-Secondary 
Y-S eccndary 



La nd 

Distance 

(Ks) 

540.730 
369. 248 
525.659 



S3 a 

Di st ance 

(lil) 

35.548 

126 .096 
15 .832 



TASC 55 
Total 
Dist ance 

(!sa) 

576. 278 
495.344 
541.491 
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An attempt was made to determine if this discrepency 
was due to program LOPLC, A comparison of rhe total 
distances in the report between the transmitters and TASC 55 
were compared to the total distances computed by INVEE1. 
The difference in distance between Subroutine INVER1 and the 
total distances from the Semi-Empirical TD 3rid acticle 
ranged from 160 - 200 m corresponding to 0.5 to 0.65 usee. 

This results in time difference errors of 0.04 usee for the 
9940-X rate and a 0. 10 ;isec for the 9940-Y rate, implying 
that the offsets, 0.86 psec and 1.13 ^sec, are caused in 
part by the method in which the azimuth from north is derer- 
mined. Unfortunately, azimuth data from TASC 55 was not 
presented in the article. The Ar^alytic Science Corporation 
has been contacted on numerous occasions in an attempt to 
ascertain their method of determining distance and azimuth. 
As of this date there has been no response. 

Nevertheless, the data from Honterey Bay was 
utilized in Program LOPLC to obtain results that could be 
compared to that obtained by the other methods. If this 

method is accurate enough, ASF Correctors could be deter- 
mined via computer for each individual position without 
using tables or field dersrmined correctors. The mean 
offset (x) and standard error (s) in microseconds between 
the observed and calculated rates were computed using equa- 
tions 3.1 and 3.2, respectively. The mean offset (1) and 

drms in meters were computed using equation 3.3 and 3.4, 
respectively. 

3 • Calculated Table ASF Corr ectors 

The offset between the observed LORAN-C rates and 
the expected values with applied AS? Correctors from the 
LORAN-C Correction Tables and the seawater Secondary Factor 
(equation 2.13) were determined using program LORTAB which 
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was written by the author (Appendix D) . The ASF Corrector 
for each position was selected using subroutine TABLE. 
Subroutine TABLE, which is found in program LORTAB deter- 
mines an ASF Corrector for each data point based on its 
geodetic position. The ASF Correctors used in subroutine 
TABLE (see Figures 2.5, 2.6, and 2.7) are located between 

latitude 36o 35' N and 36° 55’N and longitude 121° 50* W and 
122° 00’ W. The difference in sign between the ASF 

Correctors in subroutine TABLE and those found in the 
LORAN-C Correction Table is due to the difference in their 
application. AS? Correctors from the tables are applied to 
observed rates while ASF Correctors from Subroutine TABLE 
are applied to the calculated LORAN-C rates. Negative ASF 
Correctors from the LORAN-C Correction Table were applied to 
the calculated time differences to be consisrent with the 
application of Secondary Factors to the computed primary 
phase delay. 

As before equations 3.1 through 3.4 were used to 
compute -he the mean offset , standard error, and drms. If 
this application of LORAN-C Correction Tables is accurate 
enough, it precludes the need to determine ASF Correctors in 
the field. 

• Observed ASF Corre c tors 

Observed ASF Correctors were determined using 
Program ASFSEL (ASF Selection) which was written by the 
author (Appendix E) , Schnebele's prior data [1979], and the 
June 1982 data. Program ASFSEL (ASF Selection) was written 
by the author. This program computes the ASF Correctors by 
subtracting the observed LORAN-C rates from the expected 
values. Only the seawater Secondary Factors from equation 
2.13 have been applied to the calculated time differences. 
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The mean ASF Correctors for the LORAN-C rates were deter- 
mined at every minute of latitude and longitude between 36^ 
50' N and 360 35* N and 122° 04*ff and 122° 49* M. See 
Appendix E for mean ASF Correctors at one minute intervals. 
The ASF Correctors were then selected and assigned to 
subroutine TABLE in Program LORTAB at five minute latitude 
and longitude intervals. All values were enrered to the 
nearest hundredth of a microsecond. Equations 3.1 through 
3.4 were used as before for computations. This determina- 
tion cf ASF Correctors in the field, if accurate enough, may 
allow the use of LORAN-C as a positioning system for hydro- 
graphic surveying in the future. 
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IV 



RESULTS 



A total of 620 time differences (TD) and geographic 

positions based on four lines of position were recorded in 
the southern portion of Hontery Bay in order to compare the 
use of differential LORAN-C with ASF Multiple Correctors. 
The data was divided into four sets. The first data set is 
Schnebele's [1979] prior data consisting of a total of 130 
data points collected on two separate days, June 12 and July 

25, 1979 between latitude 360 38* N and 36° 47* N and 

longitude 121® 49* W and 122° 02* W. The recorded LORAN-C 

time difference rates were 9940-Y and 9940-W. 

The second data set consists of 193 time differences 
collected on June 7, 1982. The recorded LORAN-C rates were 

also 9940-Y and 9940-W. Data set Number 2 is located 

between latitude 36o 40* N and 36<> 45* N and longitude 121° 
54’ W and 122° 00* W. See Figure 3.1 for the location of 
test areas in Monterey Bay, California. 

Data set Number 3 is located in the same area as data 
set Number 2, between latitude 36o 40* N and 36o 45* N and 
longitude 121° H 54* and 12 20 00* W. This set, consisting 
of 128 points with recorded rates 9940-X and 9940-Y, was 
recorded on June 6, 1982. 

The final set, data set Number 4 was recorded between 
June 3 and ,June 5, 1982. It contains 169 points located 

between latitude 36° 36 H and 36° 39* N and longitude 121® 
53* W and 121° 58* W. The recorded rates were 9940 -X and 
9940-Y which are the same as data set three. The dara for 
rates 9940-X and 9940-Y was Icepr separate so as to determine 
if there was a significant difference between the offshore 
(data set three) and inshore (data set four) drms values due 
to phase recovery (see Chapter Two) . 
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TABLE X 



Dara Set Parameter 



Set 

Number 


Collection 

Date 


TD 

R ates 


Number of 
Sail Points 


Area 


Limits 

Lon 


1 


June 

July 


12, 1979 
25, 1979 


9 940- Y 
9 94 0- W 


130 


36/38 

36/47 


121/49 

122/02 


2 


June 


7, 1982 


9 940- Y 
9 940-W 


193 


36/40 

36/45 


121/54 

122/00 


3 


June 


6, 1982 


9 94 0- X 
9 940- Y 


1 28 


36/40 

36/45 


121/54 

122/00 


4 


June 


5, 1982 


9 940- X 
9 940-Y 


169 


36/36 

36/39 


121/53 

121/58 



i I 



Table X provides a convenient breakdown of the parame- 
ters for each of the data sets. The table consists of the 
data set number, the date the data was collected, the 
LORAN-C time difference rates, the number of data poinxs, 
and the area limits in latitude and longitude. It defines 
the parameters for the four data sets of Tables XI (Seawarer 
Secondary Factor Errors) , XII (Semi-Empirical TD Grid 
Correction Errors), XIII (Multiple LORAN-C Correc-ion Table 
Errors) , and XIV (Multiple Observed Correction Errors) . All 
basic data and computarions are kept on file with the NOAA 
hydrography instructor at the Naval Postgraduate School, 
Monterey, California. (Individuals seeking this information 
should contact the Oceanography Department.) 

A. SEAWATER SECONDARY FACTORS 

LORAN-C rates were computed using only the seawater 
Secondary Factor (equation 2.13). The LORAN-C computed time 
differences are basically uncorrected rates since no 
Additional Secondary Factors were applied. Offsets between 
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TABLE XI 



Seawater Secondary Factor Errors 









Mean 


Mean 


Standard 










Offset 


Offset 


Error 


drms 


Rates 




fx-jisec) 


azi) 




(1) 


4) 


Set 


1 


(Schnebele ' s 


da t a , 130 


data points) 






Y 




-0.508 


-1 48.636 


t).088 


110.4 




W 




- 1. 24 1 


-6 83.693 


0.134 




b) 


Set 


2 


(June 7, 1982 


, 1 93 data 


points) 






Y 




-0.526 


-1 54.295 


0.077 


101.3 




W 




-1.283 


-701.947 


0.124 




c) 


Set 


3 


(June 6, 1982 


, 128 data 


points) 






X 




1.565 


5 35.449 


0.059 


38.6 




Y 




-0.550 


-1 61.404 


0.083 




d) 


Set 


4 


(June 3-5. 1982, 169 da 


ta points) 






X 




1.582 


5 50.089 


0.055 


68.1 




Y 




-0.817 


-237.109 


0. 187 





the observed and calculated rates using only the seawater 
Secondary Factor were computed to illustrate rhe improvement 
in position after applying the Semi-Empirical TD Model, ASF 
Loran-C Correction Tables, or the Multiple Observed ASF 
Correctors. The mean offset in microseconds and meters, 
standard error, and the drms are found in Table XI. 
Examples of offsets for several data points are listed in 
Appendix 3 after program LORAN. 

It was stated earlier in Chapter One that Schnebele 
obtained a 66 m 1 drms using Y and W rates. The drms of 66 
m was obtained using 48 data points which were located 10 km 
or more offshore. The 110.4 m 1 drms for the Y and W rates 
in Table XI is a result of combining the 130 inshore and 
offshore positions. The increase from 66.0 m 1 drms for the 
offshore positions to 110.4 m for combined offshore and 
inshore positions indicates that the application of a single 
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ASF Corrector is dependant on the size of the area. The 48 
data points were located between latitude 36° 41 » N and 36° 
46' N and longitude 121° 55 » W and 122° 02’ W, 

Schnebele also obtained the 66 m 1 drms by adjusting the 
LORAN-C observed time differences which were skewed due to 
the ship's motion and the five second averaging interval of 
the LORAN-C receiver. These caused the observed TDs to be 
several seconds old in comparison to the microwave system 
measurements [Schnebele, 1979]. Due to the large amounrs of 
data from the June 1982 survey operations, no deskewing was 
done. 

The offset, standard error, and drms for data set Number 
2 (June 7, 1982) compares well to data set Number 1 

(Schnebele's data - all). Also, the 38.6 m 1 drms for the 
X-Y rates is between Schnebele's 42.0 m 1 drms prediction 
[Schnebele, 1979 ] and Nelson's findings of 30 ra 1 drms in 
San Francisco Bay [General Electric Go., 1979]. The large 
drms of 68.1 m for data set Number 4 (X-Y rates) is probably 
due to the phase recovery of the electromagnetic wave from 
The 9940-Y transmitter located in Searchlight, Nevada. 
Finally, if ASF CorrecTors are not applied, drms values 
ranged from 38.6 m for 9940 X-Y to 101.4 m for 9940 Y-W for 
data sets 2 and 3. The large offsets for all data sets 
indicate a systematic error, ranging from 150 m to 700 m, 
which precludes the use of this method for hydrographic 
surveying. 

3. SEMI -EMPIRICAL TD GRID 

To determine if the Semi-Empirical TD Grid would reduce the 
offset and drms for hydrographic surveying, program LOPLC 
(Appendix C) was applied to dara sets 1 through 4. The 
Semi- Empirical TD Grid applies a Secondary FacTor and an 
Additional Secondary Factor to the primary phase delay based 



57 



and the 



on the distance overland, the distance over water, 
total distance using mean land and water conductivities. 
The mean offset in meters and microseconds, standard error, 
and drms are listed in Table XII. Examples of offsets for 
several data points can be found in appendix C after program 

TABLE XII 



Semi-Empirical ID Grid Correction Errors 



Rates 



Mean 

Offset 

(xr^sec) 



Mean 

Offset 

(iia) 



standard 

Error 

( s-pse c) 



drms 

(a) 



a) 


Set 


1 


(Schnebele ' s 


dat 


a, 1 


30 : 


iat 


a 


points) 






y 




1.131 




330. 


672 






0. 


086 


113. 


,5 




w 




1.145 




630. 


278 






0. 


139 






b) 


Set 


2 


(June 7, 1982 


da 


ta, 


193 


da 


ta 


po 


ints) 








Y 




1. 108 




324. 


975 






0. 


077 


112, 


,7 




w 




1. 105 




603. 


926 






0. 


126 






c) 


Set 


3 


(June 6, 1982 


da 


ta. 


128 


da 


ta 


po 


ints) 








X 




1.131 




38^ 


940 






0. 


061 


39. 


.4 




Y 




0.5 89 




172. 


923 






0. 


084 






d) 


Set 


4 


(June 3-5, 19 


82 


data 


, 169 


da 


ta 


points) 






X 




1.0 5iD 




365. 


113 






0. 


052 


63. 


.2 




Y 




0.341 




99. 


084 






0. 


189 







LOPLC. Earlier, Program LOPLC had been tested with data 
found in the article by The Analytic Science Corporation 
[1979]. Program LOPLC produced offsets of 0.86 ps^c and 
1.13 psec for the time differences from rates 9940-X and 
9940-Y, respectively. As previously stated, the large 
offsets may be due tc the method by which the distance and 
azimuth were computed. The drms, for the data from Monterey 
Bay, obtained with the Semi -Empirical TD Grid was similar to 
the drms errors for seawater Secondary Factor Model. This 
might be an indication that Program LOPLC is correct but the 
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bias needs to be adjusted to reduce the large offset between 
the observed and calculated ID rates. The mean offset 
ranged from 99 to 630 m. Bias reduction could be achieved by 
applying land-sea data for the entire West Coast to rhe 
Semi-Empirical model. Again, the existence of large offset 
values precludes the use of this method for hydrographic 
surveying. 

C. TABLE ASF CORRECTOHS 

To determine if the multiple ASF corrections from the 
LORAN-C Correction Table would reduce the offset and drms to 
meex the NOS accuracy standards, program LORTAB (Appendix D) 
was applied to data sets 1 through 4. ASF Correctors from 

TABLE XIII 

Multiple LORAN-C Correction Table Errors 



( 1 









Mean 


Mean 


standard 










Offset 


Offset 


Error 


drms 




Rates 


(x-psec) 


(iza) 


(si;^sec) 


(S) 


a) 


Set 


1 


(Schnebele ' s 


data, 130 ; 
-67. 068 


iata points) 


123.6 


y 




-0.229 


0. 104 




w 




0.169 


92. 540 


0. 160 




b) 


Set 


2 


(June 7, 1982 
-0.309 


data, 192 
-90. 614 


data points) 




Y 




0.085 


116.9 




W 




0.123 


67. 189 


0. 144 




c) 


Set 


3 


(June 6, 1982 
0.559 


data, 128 
191. 134 


data points) 




X 




0.073 


47.3 




y 




-0.319 


-93. 486 


0. 101 




d) 


Set 


4 


(June 3-5, 1982 data, 169 data points) 




No 


AS 


F Corrections 


listed in 


Table for the 


south- 



east end of Monterey Bay next zo rhe shore line. 



the LORAN-C Correction Tables are determined from field 
observation of land conductivities by the U.S. Coast Guard 
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Calibration Team [Marine Science Department, 
results are listed in Table XIII. 



1982]. The 



When compared to the Seawater Secondary Factor Error in 
Table XI, rhe drms values using LORAN-C Correction Tables 
were increased slightly while the offsets were reduced 
substantially. The offsets ranged from 65 - 200 m. Since 
the drms ranged from 47.3 to 116.9 m for rates 9940 X-Y and 
9940 Y-W respectively, the application of ASF Corrector from 
the tables does not meet the NOS accuracy standard. Again, 
large offsets and the increase in drms precludes the use of 
this method for hydrographic surveying. 

D. MULTIPLE OBSERVED CORRECTORS 

To determine if multiple observed correctors would 
diminish the offset and drms values, mean ASF Correctors 
were selected at one minute latitude and longitude intervals 
using Program ASFSEL. One minute ASF Correctors are shown 
at the end of Program ASFSEL in Appendix E. From the one 
minute grid, mean ASF Correctors were selected and entered 
into Subroutine TABLE at five minute latitude and longitude 
intervals in Program LORTA3 . The following is an example of 
the 9940-X ASF Correctors at five minute latitude and longi- 



tude intervals 


for the prog 


ram: 






122/00/00.0 


1 21/55/00.0 


121/50/00.0 


36/50/00.0 


-1. 52 


-1.63 




36/45/00.0 


-1. 52 


-1.56 


-1.61 


36/40/00.0 




-1.58 


- 1.60 



This is the same format used in the LORAN-C Correction 
Tables. Program LORTAB was applied to the four data sets. 
The error results are listed in Table XIV. 

The drms value obtained with multiple observed correc- 
tors were all reduced when compared to the arms for nhe 
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7^3 




rABLE XIV 



Multiple Observed Correction Errors 



Rates 



a) 



b) 



c) 



d) 



Mean 
Offset 
( X-fl sec ) 



Mean 

Offset 

1-m) 



Standard 

Error 

(s-jasec) 



Set 


1 


(Schnebel! 


g 


' s da 


t 


a, 130 


dat 


a 


point: 


Y 






0. 


033 








9.585 






0.076 


W 






0. 


052 






2 


8. 604 






0.113 


Set 


2 


(d 


une 


7. 


1 


982, 


1 


93 data 


po 


in 


ts) 


Y 




0. 


028 








8.31 8 




0.073 


w 






-0. 


016 






— 


8. 944 






0.1 16 


Set 


3 


(J 


une 




1 


982, 


1 


28 data 


po 


in 


ts) 


X 




-0. 


034 






- 


11.490 




0.055 


Y 






0. 


010 








2.9 13 






0.086 


Set 


4 


(J 


une 


3.5 


/ 


1982 




169 data 


po 


ints) 


X 




0. 


001 








0.383 




0.052 


Y 






-0. 


008 








-2.225 






0.187 



drms 

(B) 



87.6 

89.3 

38.3 
67.5 



seawater Secondary Corrector. The most impressive reduction 
in drms was within Schnebele's data which covered an area of 
seven minutes of latitude and 12 minutes of longirude. The 
drms for seawater Secondary Correctors was 110.4 m whereas 
the drms for the same data using multiple observed correc- 
tors was 87.6 m. This is a smaller drms than that of the 
June 7, 1982 data (data set Number 2) which was obtained 
three years later. It appears that 87.6 m 1 drms is nearly 
the minimum error that can be obtained for the 9940 Y-W 
rates in Monterey Bay after applying multiple observed AS? 
Correctors at five minute latitude and longitude intervals. 
For the LORAN-C rates 9940 X-Y, a 38.3 m 1 drms was obtained 
for same five minute area covered by data set Number 2 (June 
7, 1982) . LORAN-C rates 9940 X-Y were not obtained for the 
same size area covered by data set Number 1 (seven minutes 
of latitude and twelve minutes of longitude - 9940 Y-W) due 
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to the restriction of ship time and the length oi time the 
Racal-DECCA Trisponder electronic eqaipment had been loaned. 

The drms for the June 3-5, 1982 inshore data was only 
reduced to 67.5 m from 68.1 m for seawater Secondary 
Correctors. The small change in error at the coast is 
probably a result of the erratic behavior of phase recovery 
discussed earlier in Chapter Two. The drnis value could 
presumably be reduced if the correctors were applied at one 
minute intervals. This would be a very costly method of 
calibrating Loran-C for hydrographic surveying. 
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V. CONCLUSIONS 

It was noted in Chapter One rhat the smallest scale 
routinely used for coastal surveys is 1 ;80,000, This yields 
an allowable error of 40 m 1 drms with no systematic errors. 
This paper determined whether or not applying multiple 
Additional Secondary Factors (ASF) Correctors to LORAN-C 
lines of position would reduce the drms sufficiently to meet 
the accuracy standards set by the National Ocean Survey. 

Three methods of applying multiple ASF Correctors were 
tested. The first approach computes the time difference 
based on a Se mi-Smpir ical TD Grid. The Semi-Empirical Model 
produced large offsets in the 9940-W, 9940-X, and 9940-Y 

time differences. The offsets ranged from 99 to 630 m. The 
drms for 9940 X-Y combination was 39.4 m and the drms for 
9940 W-Y combination was 102.7 m. 

The second method applies ASF Correctors found in the 
DMAHTC LORAN-C Correction Tables to LORAN-C lines of 
position. The application of the tables reduced the offset 
in the LCRAN-C time differences. The offsets were between 
67 and 191 m. The drms was increased to 47.3 m for the 9940 
X-Y combination and 116.9 m for the 9940 W-Y pair. 

The final and most accurate method applies multiple 
observed ASF Correctors at five minute latitude and 
longitude intervals to LORAN-C lines of position. This 



method again 


reduced the 


offset in 


the 


time difference. 


This offset 


was between 


3 and 12 


m 


for the 9940 X-Y 


com bination . 


Part of the 


offset may 


hav 


e been a result of 


the microwave 


! positioning 


system. 


R ef e 


rence is made to 



Anderson's [ 1982 ] paper (in preparation) on the evaluation of 
multiple lines of position. 
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The drms values were also reduced to 38.3 m for the 9940 
X-y rates and 89.3 m for the 9940 W-Y combination. The 38.3 
m 1 drms can be decreased by improving the sampling time for 
LORAN-C receivers. Nelson obtained 30 m 1 drms for the 9940 
X-Y rares with special LORAN-C eguipmenr used in San 
Francisco Bay, California [General Electric Co., 1979]. 
Improving rhe sampling time for LORAN-C receivers used as 
positioning equipment for hydrographic surveys should be 
investigated. 

With drms values of 38.3 m with the possibility of 
obtaining 30.0 m 1 drms and offsets ranging from 3 to 12 m, 
it may be possible to use LORAN-C for hydrographic surveys 
at scales of 1:80,000 or less using multiple observed ASF 
Correctors. The use of DMAHTC LORAN-C Correction Tables 
should not be ignored. After updating these ASF Correctors 
with observed data, the LORAN-C Correction Tables may allow 
LORAN-C to be used as a positioning system for hydrographic 
surveys. 
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APP ENDI X A 



PRO G HAM GPB YLQ 



C PROGRAM GPBYLQ 

C 

C GENERAL PROGRAM FOR DETERMINING GP FROM KNOWN STATION 

C POSITIONS AND THE DISTANCES FROM THEM DSING LEAST 

C SQUARES. GP TO XY AND XY TO GP ARE DETERMINED BY 

C SUBROUTINES GPTOXY AND XYTOGP WHICH ARE BASED ON THE 

C MODIFIED TRANSVERSE MERCATOR PROJECTION (MTM) . 

C 

C PROGRAMMED BY GERALD E. WHEATON, LT. NOAA 

C 

C LEAST SQUARES ADJUSTMENT BY PAUL R. WOLF, PH.D. AND 

C REVISED BY D, LEATH, LCDR 

C 

C PROGRAM INPUT VARIABLE NAMES 

C TITLE = ANY JOB IDENTIFICATION NAMES OR NUMBERS 

C M AND N = THE NUMBER OF EQUATIONS (M) AND UNKNOWNS (N) 

C XO AND YO = BEST QUES TIMATE OF THE POSTION 

C STA(I,1) = X COORDINATE OF KNOW STATION 

C STA(I,2) = Y COORDINATE OF KNOW STATION 

C STA(I,3) = DISTANCE FROM KNOW STATION 

C A(I,J) = THE COEFFICIENT MATRIX 

C EL(I,J) = THE CONSTANT MATRIX 

C QLL(I,J) = THE WEIGHT MATRIX (WEIGHTS ARE ENTERED AS 

C 1'S IF THE SOLUTION IS EQUALLY WEIGHTED) 

C 

IMPLICIT REAL+8 (A-H,0-S) 

DIMENSION XCORD (30) ,YCORD (30) ,WT (30) ,ISNO(30) 

COMMON /ISTAT/ STA (30 , 3) , FXY (30) 

COMMON /L3QX/ A (30 ,30 ) , EL ( 30 , 1) , QLL ( 30 , 30) , AT ( 30 , 30 ) , 
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1AQ(30,30), QXX(30,30) , AQL(30,10), X(30,1), V{30,1), 
2VAH(30) , TITLE(80) 

C 

C READ AND WRITE OUTPUT TITLE 

C 

WRITS (6,509) 

HEAD (5,71 0) TITLE 
WRITE (6,710) TITLE 
C 

C DEFINE NUMBER OF SI3NALS (NOT GREATER THAN 30) AND 

C NUMBER OF DATA SETS. 

C 

READ (5, 502) NSIG,NDATA 

502 FORMAT(I5,I10) 

C 

C DEFINE THE CENTRAL MERIDIAN (CMER) , FALSE EASTING 

C (FEST), AND CENTRAL LATITUDE (CLAT) . 

C CMER IS EXPRESSED IN DEGREES, MINUTES, AND SECONDS. 

C FEST IS THE X-COORDINATE THAT IS ASSIGNED TO THE 

C CENTRAL MERIDIAN AND IS EXPRESSED IN METERS. 

C CLAT IS DEFINED AS THE CONTROLLING LATITUDE. 

C IT IS USED TO REFERENCE THE Y-COORDIN ATES AND 

C IS EXPRESSED IN METERS. 

C 

READ (5,503) ILO NC, ILM INC, RLSECC , FEST ,CLAT 

503 FORMAT(1X,I3, 1X,I2,1X ,?8. 5, 1X,F7. 1 ,?10. 1) 

CMER = ( (IA3S (ILONC) * 60 + ILMINC) * 60) + RLSECC 

C 

C DEFINE VARIABLE FOR; 

C NUMBER OF EQUATIONS (M) 

C NUMBER OF UNKNOWNS (N) 

C IPAGE = NUMBER OF LINES PER PAGE. 

C 



IPAGE 



1 
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c 



H = 4 
N = 2 



C NOLL WEIGHTS 

C 

DO 3 I M= 1,30,1 
DO 3 JM=1,30,1 

3 QLL(IM,JM) = 0.0 

C 

C READ STATION NUMBERS , POSITION, AND WEIGHTS. 

C CONVERT POSITIONS (3P) TO SECONDS AND THEN TO XY. 

C 

DO 12 J=1,NSIG,1 

READ (5,800)ISNO (J) ,IL AT ,I MIN , ES EC, JLON, JHIN , SS EC , WT (J) 
800 FORM AT (IX, 1 3, 14, 13, F9 .5 ,1 5 , 13 , F9 . 5 , F5 . 1 ) 

C 

RMAST = ((lABS(ILAT) ♦ 60 + IMIN) * 60) + KSEC 
RMASTL = ((lABS(JLON) * 60 + JMIN) * 60) + SSEC 

C 



12 

C 

C 

c 

15 



CALL GPTOXY (RMAST, RMASTL, XMETER , YMET ER, FEST , CLAT , CMER) 
XCORD(J) = XMETER 
YCORD(J) = YMETER 

READ DATA (STATION NUMBERS AND THE DISTANCES) 

DO 40 JCOUNT = 1 ,NDATA , 1 



805 


READ(5, 805) IF, 
FORM AT (515) 

DO 16 ICODNT=1 


IS, IT, 
,NSIG, 


14, 

1 


NREC 




IF (ISNO (ICOUNT) 


.EQ. 


IF) 


IF=ICOUNT 




IF (ISNO (ICOUNT) 


,EQ. 


IS) 


IS=ICOUNT 




IF (ISNO (ICOUNT) 


.EQ. 


IT) 


IT=ICOUNT 


16 

C 


IF (ISNO (ICOUNT) 


.EQ. 


14) 


I4=IC0UNT 
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c 



c 



c 



c 

c 

c 



507 

C 

C 

C 

c 

19 
C 

c 

c 

c 

20 



STA(1,1) = 
STA(1,2) = 
QLL(1,1) = 

STA(2,1) = 
STA(2,2) = 
QLL(2,2) = 

STA(3,1) = 
STA(3,2) = 
QLL(3,3) = 

STA(4,1) = 
STA(4,2) = 

QLL(4,4) = 



XCORD (I?) 
YCOED (IF) 
WT (IF) 

XCORD (IS) 
YCORD(IS) 
WT (IS) 

XCORD (IT) 
YCORD(IT) 
MT (IT) 

XCORD (14) 
YCORD(I4) 
WT (14) 



READ THE DISTANCE RECORD AND LORAN HATS. 



DO 38 KCOONT=1 , NREC, 1 

READ (5,507) STA (1, 3) , STA(2,3) ,STA (3,3) ,STA (4,3) , 
1 EATS1,RATE2 
FORMAT (6F10. 1) 

DETERMINE BEST GUESS COORDINATES XO AND YO 
WITH SUBROUTINE GUESS. 



CALL GUESS (XO,YO) 
IJUM? = 0 



COMPUTE FXY, A AND L MATRIX 
DO 25 1=1, M 

FXY(I) = DSQRT (DABS ( ( XO-STA (I, 1) ) *<=2 
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1 (YO-STA (I, 2) )**2) ) 

A (1,1) = (XO-STA(I,1) ) / (FXY(D) 

A(I,2) = (YO-STA(I,2) ) / (FXY (I) ) 

25 SL(I,1) = STA(I,3) - FXY(I) 

C 

C CALL SUBROUTINE LSQR 

C 

CALL LSQR(M,N) 

C 

C COMPUTE THE NEW QUESSTIMATE FOR XO AND YO 

C 

XO = XO + X (1 , 1 ) 

YO = YO + X (2,1 ) 

C 

C EXIT IF STANDARD ARE MET USING IJUMP OR 

C XO AND YO CUT OFF 

C 

IF (DABS (X (1 , 1) ) .LE. 1 . 00 .AND. DABS (X (2, 1) ) . LE. 1 . 00) 
1 GO TO 35 

IF (IJUMP .EQ. 10) GO TO 3 5 
30 GO TO 20 

C 

C COMPUTE ERROR ELLIPSE 

C 

35 CONTINUE 

CALL ELIPSE (SO, SV) 

C 

C CONVERT XY TO GP 

C 

CALL XYTCGP (XO, YO, SEC LAT, SECLON , F SST , CLAT , CM ER) 

CALL TODMS (SECL AT, IDE GP ,I MI HP, RS ECP) 

CALL TODMS (SECL ON, JDE GP , J MINP , 3 SEC P) 

C 

C PAGE AND CONTINUE WITH NEXT SET OF OBSERVATIONS 
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c 

WRITE (6,5 05) ID EGP ,IM INP,RSECP, JDEGP, JMINP , SSECP , 

1 RATE1,RATE2 

C WRITE (6,505) ID EGP ,IH INP, RSECP, JDEGP, JMINP , SSEC? , 

1 SO,SV 

505 FORMAT(I4,I3,F6.2,I5, 13 ,F 6 . 2 , 2F9 . 2) 

C IF(IPAGE .EQ, 50) WRITE (6,509) 

IF(IPAGE .EQ. 50) IPAGE = 0 
38 IPAGE = IPAGE + 1 

C 

C FORMAT STATEMENTS 

C 

710 FORMAT(80A1) 

509 FORMAT(IHI) 

40 CONTINUE 

STOP 
END 

SUBROUTINE LSQR(M,N) 

C 

IMPLICIT REAL*8 (A-H,0-Z) 

COMMON /LSQX/ A (30 , 30 ) , EL (30, 1 ) , QLL (30 , 30) , AT ( 30 , 30) , 
1AQ(30,30), QXX(30,30) , AQL(30,10), X(30,1), V(30,1), 
2VAR(30) , TITLE (80) 

C 

C COMPUTE A TRANSPOSE 31 TRANSPOSING THE A MATRIX (AT) 

C 

DO 61 1=1, M 
DO 6 1 J=1 ,N 
61 AT (J ,1) =A (I,J) 

C 

C USING STEPS (1) , (2) , AND (3) COMPUTE THE INVERSE 

C OF THE TRANSPOSE (AT) ♦ WEIGHTED MATRIC (QLL) * 

C MATRIX A = QXX. 
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c 

C (1) COMPOTE AQ = AT * QLL 

C 

DO 7 1 1 = 1 ,N 
DO 7 1 J=1 ,M 
AQ (I, J) =0. 

DO 7 1 K = 1 ,M 

71 A<2 (I, J) =AQ (I, J) +(AT(I ,K)*QLL(K, J) ) 

C 

C (2) COMPUTE QXX = AQ ♦ A 

C 

DO 01 1 = 1, N 
DO 81 J=1,N 
QXX(I,J)=0. 

DO 81 K=1,M 

81 QXX(I,J)=QXX (I, J) +AQ (I,K) *A (K, J) 

C 

C (3) INVERT QXX MATRIX 

C 

DO 307 K=1,N 
DO 302 J=1,N 
IF (J-K) 304,302, 304 

304 QXX(K,J) =QXX (K, J) /QXX (K,K) 

302 CONTINUE 

QXX(K,K) = 1./QXX (K,K) 

DO 307 1=1, N 
IF (I-K) 305,307, 305 

305 DO 303 J=1,N 

IF (J-K) 306,303, 306 

3 06 QXX(I,J)=QXX (I, o) - QXX (I,K) *QXX (X, J) 

303 CONTINUE 

QXX(I,X)=-QXX (I,K) *QXX (K, K) 

307 CONTINUE 
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C USING STEPS (4) AND (5), COMPUTE THE UNKNOWNS X 

C BY MULT THE INVERSE Q XX AND AQL. 

C 

C (4) COMPUTE AQL = AQ * EL 

C 

DO 101 1=1, N 
AQL(I,1)=0. 

DO 101 K=1,M 

101 AQL(I,1)=AQL (I, 1) +AQ (I,K) *EL(K, 1) 

C 

C (5) COMPUTE X = QXX =»■ AQL 

DO 201 1=1, N 
C 

X (I, 1) =0. 

DO 201 K=1,N 

201 X (I, 1) =X (I, 1) +QXX (I,K) *AQL (K, 1) 

C 

C (6) COMPUTE THE RESIDUAL {V = A * X -EL) 

C 

DO 3 01 1=1, M 
V (I, 1)=0, 

DO 301 K=1,N 

301 V (I, 1)=7(I, 1) +A (I,K) *X(K, 1) 

DO 1 1=1, M 

1 V (I, 1)=V(I,1) -EL(I,1) 

C 

C COMPUTE THE STANDARD DEVIATION OF UNIT WEIGHT SIGMA 

C DM - NUMBER OF OBSERVATIONS 

C DN - NUMBER OF KNKNOWNS 

C 

SIGHA=0. 

DM=M 

DN=N 

DO 382 1=1, M 
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382 SIGMA=SIGMA+V (1,1) **2*QLL(I,I) 

SIGMA=DSQBT (SIGMA/ (D M-DM) ) 

C 

C COMPUTE THE STANDARD DEVIATION OP THE ADJUSTED UNKNOWNS 

C QXX - ARE THE ELEMENTS OF THE COVARIANCE MATRIX. 

C 

DO 446 1=1, N 

446 VAR(I) =DSQRT (QXX (I ,1) *SIGMA**2) 

C 

510 CONTINUE 
RETURN 
END 

C== ======= =========== ===^== ========== = ========== = ============ 

SUBROUTINE ELIP SE ( SU , SV) 

C 

C 

C SOLVE FOR THE SEMI MAJOR AND SEMI MINOR AXIS OF 

C THE ERROR ELLIPSE 

C 

IMPLICIT REAL*8 (A-H,0-Z) 

COMMON /LSQX/ A (30 ,3 0 ) , EL ( 30, 1 ) ,QLL (30,30) ,AT(30,30) , 
1AQ(30,30), QXX(30,30) , AQL(30,10), X(30,1), V(30,1), 
2VAH(30) , TITLE(80) 

C 

SUS = .5*(QXX(1,1) + QXX(2,2) + DSQRT (DABS (QXX ( 1 , 1) - 

1 QXX (2,2) +4.0*QXX (1,2) *QXX (2, 1) ) ) ) 

SVS = .5*(QXX(1,1) + QXX(2,2) - DSQRT (DABS (QXX (1 , 1) - 

1 QXX (2,2) +4.0*QXX (1,2 ) *QXX (2, 1) ) ) ) 

C 

SO = DSQRT (SOS) 

SV = DSQRT (SVS) 

C 

RETURN 

END 



73 



c==================== ====== =============== === = == = = == = ====== 

S03R0UTINE GUESS(XO,YO) 

C 

C 

C SUBROUTINE GUESS DETERMINES THE BEST GUESS COORDINATES 

C TO BE USED IN SUBROUTINE LSQR. USE RIGHT SIDE RULE FOR 

C STATION ORDER. 

C 

IMPLICIT REAL*8 (A-H,0-Z) 

CCHHON /ISTAT/ STA {30 , 3) , FXY (30) 

C 

C DETERMINE DISTANCE BETWEEN STATION 1 AND STATION 2 

C 

D=DSQRT ((STA (2, 1) - ST A (1 , 1 ) ) **2+ (ST A (2 , 2) -ST A (1 , 2) ) ♦ *2) 

C 

C DETERMINE ANGLE ALPHA BETWEEN XO , YO/STA2/ST A 1 

C 

ALPHA=DARCOS ( (STA (2,3 ) ♦♦2- STA (1,3) «*2 + D**2) / 

1 (2. 0*STA (2,3) *D) ) 

C 

C DETERMINE ANGLE BROVO BETWEEN X-AXIS AND STA2-STA1 

C 

BROVO=DARSIN ( (STA (1,2 ) -ST A (2,2) ) /D) 

C 

C DETERMINE X AND Y LENGTH 

C 

X=STA (2,3) ♦DCOS (ALPHA +BRO VO) 

Y=STA (2,3) *DSIN (ALPHA +BROVO) 

C 

C DETERMINE XO AND YO 

C 

XO = STA(2, 1) +X 
YO = STA (2,2) +Y 
C 
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HETURN 

END 

C== ============ ====== ====== ======= === = ==== ====== = === = ======== 

SUBROUTINE GPTO XY ( SEC LAT, S ECLON , XCO, YCO, FEST , CL AT ,CMER) 

C 

C 

IMPLICIT REAL*8 (A-H,0-Z) 

DATA E2^RKO,A / .00 676 8658 DO 99 99 8 DO , 63782 0 6 . 4 DO/ 

DATA RKGEO,W1,W2 /O. 0 4848 1 368D0 , 0 . 1 1 4 22D0 , 2 1 . 7360 7D 0/ 
DATA W3,W 4/51 04 .57 33 8 DO, 63 673 99.6 39D0/ 

DATA RADSEC /. 0 000 04 3 48 13 68 1 1 1 DO / 

C 

C 

RADLAT = SECLAT * RADSEC 
SINLAT = DSIN (RADLAT) 

SIN2LA = SINLAT * SINLAT 
COSLAT = DCOS (RADLAT) 

COS2LA = COSLAT ^ COS LAT 
P = (CMER - SECLON) / 10000. ODO 
V = A / DSQRT(1.0D0 - E2*SIN2LA) 

TANCON = 1.0D0 - SIN2 LA/COS2LA 
S = W4 * (RADLAT - SI NLAT*COSLAT/1 0. 0D0**6 * 

* (W3 -COS2LA* (W2- W1*COS2LA) ) ) 

T1 = S * RKO - CLAT 
T2 = RKGEO ♦ CO SLAT * RKO * V 
T3 = T2 * RKGEO / 2. ODO 
T4 = T3 * SINLAT 

T5 = T3 * RKGEO * COS 2LA / 3. ODO 
T6 = T5 * TANCON 

T7 = (4. ODO + TANCON) *T5 *RKGEO * SINLAT / 4. ODO 

XCO = (T2 + (T6*P**2) )*P + FEST 

YCO = (T7=«‘P**4) + (T4*P#*2) + T1 

RETURN 

END 
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c== ========================================================== 

SUBROUTINE XYTOGP(XCO , YCO , SECLAT , SECLON , FE ST , CLAT , CHER) 

C 

C 

IMPLICIT REALMS (A-H,0-Z) 

C 

DATA E2^A,SR /O ,00676 8658D0 ,6378206. 4D0 , 

1 0.0000048U81368D0/ 

DATA W1,W2,W3 /O. 246 8 2D0, 30. 023 3 5D0 , 5 078 .64 977D0/ 

C 

D = CLAT + YCO 

MO = 0. 15704998 1D0/10.0D0**6 * D 
SINWO = DSIN(WO) 

COSMO = DCOS (WO) 

COS2WO = COS WO * COSWO 

PHI1 = WO + SIN WO*COS WO/1 0.0D0**6 * 

* (W3+COS2WO*{W2+W 1*COS2WO) ) 

PHI2 = PHI1 / 0.99993D0 

PHI3 = PHI2 / SR 

Q = (XCO - FSST) / 10.0D0**6 

V = A / DSQRT(1.0D0 - E2*DSIN (PHI2) **2) 

T = DCOS(PHI2) * SR 
C = V * 0.99998D0 
T1 = 10.0D0**6 / (T*C) 

T2 = (T1*10. 0D0**6) / (2.0D0*C) 

T3 = (T2*10. 0D0**6) / (3. 0D0*C) 

T4 = (T3*10. 0D0**6) / (4.0D0*C) 

DSLLON = (T1-Q=i‘*2*T3* (2.0D0’i‘DTAN (PHI2) **2+ 1 . ODO) ) *Q 
SECLAT = (( (3.0D0*DTAN(PHI2) **2 + 5.0D0) *T4) *Q=i'*2-T2) 

1 ♦Q*=«'2*DSIN (PHI2) + PHI3 

SECLON = CMEH - DELLON 
RETURN 
END 

SENTRY 
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DATA SET EXAMPLE 



PROGRAM GPBYLQ 



SHIP’S POSITION OBSERVED LORAN RATES 



LATITUDE 

(D-M-S) 


LONGITUDE 

(D-M-S) 


9940-X 

(fisec) 


9940-Y 

(^sec) 


36 


36 


40. 20 


121 


52 


48.62 


27508.79 


42742.71 


36 


36 


42.41 


121 


52 


48.16 


27508. 8 1 


42742. 93 


36 


36 


44. 56 


121 


52 


47.52 


27508.83 


42743.15 


36 


36 


46. 84 


121 


52 


46.98 


27508.76 


42743. 26 


36 


36 


49.09 


121 


52 


46.59 


27508.76 


42743.48 
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AP PBNDI X B 
PROGRAM LORAN 



C PROGRAM LORAN 

C 

C PROGRAM COMPOTES LINS OF POSITIONS FOR LORAN-C USING 

C THE SECONDARY FACTOR (SF) BASED ON SEAWATER EM MODEL. 

C 

IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION AXIS ( 13) ,RF (13) 

DIMENSION XSLAV 1 (1000 ) ,XSLA72 (100 0) 

DATA RHOSEC,PI, UNCOV, RN/2 . 06 264 80 6 24 7 1 DOS , 

1 3.1415926535898D0,299,792458D0,1 .00033800/ 

DATA XMEAN1 ,XMEAN2,VA R1,VAR2/0. OODO, 0. OODO , 

1 0.00D0,00.0D0/ 

DATA XMEAN3,XMEAN4/0. OODO,O.OODO/ 

C 

CCC CCC**** ********* *#♦*♦:)£******* 

DATA AXIS/6. 3782064006,6. 378388 D06 , 6 . 377397 1 55D0 6 , 

1 6.37816 D06, 6. 378160 06,6. 378249145D06, 6. 378165D06, 

2 6.378166006,6.378165 D06 , 6 . 378 1 4 5D0 6 , 6 . 3775634D06, 

3 6.378245006,6.3781350006/ 

DATA RF/6. 3565838 D06, 2. 97 D02 , 2 . 99 15 28 1285D02 , 

1 2.98 2500 2,2. 98 247 16 7 427D 02, 2. 93 465D 02 , 2 . 98 25 D02 , 

2 2.983 D02,2.983D02,2.9825 D02, 6. 3562569D06 , 

3 2.983 D02,2.9826D02/ 

C*****ELLIPSOID OPTION NUMBER 



c 


1. CLARKE 1866 


8. 


MERCURY 


c 


2. INTERNATIONAL (HAYFORD) 


9. 


MARSHALL ISLAND 


c 


3. BESSEL 1841 


10. 


NAVY 8D 


c 


4. AND (AUSTRALIAN) 


1 1. 


AIRY 


c 


5. 1967 REFERENCE 


12. 


KRASSOWSKI 1940 
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C 6. CLARKE 1880 MOD 13. WGS 1972 

C 7. SAO 

C 

CCC CCC** ♦♦**♦♦♦* ♦*♦♦♦♦ 4t:«c4c4( :(c4c4c:4(4c4c4c4(4c 

C CC1-2 = ELLIPSOID NUMBER (K) 

C CC3-5 = NUMBER OF POINTS ALONG COAST (IREC) 

C 

READ (5, 100) K,IREC 
100 FORMAT ( IX, 12, 14) 

TWOPI=2.*PI 
A=AXIS (K) 

F=1. /RF (K) 

IF (F.LT.3.D-3) F=(A-1./E)/A 

C 

C READ MASTER AND SLAVE STATIONS POSITIONS 

C THE FIRST RECORD IS THE NUMBER OF MASTER AND SLAVE 

C STATIONS FOR THE PARTICULAR CHAIN. 

C 

WRITE (6,202) 

202 FORMAT(IHI) 

WRITE (6, 201) 

201 FORMAT (1H ) 

READ (5, 105) ILATM,IMINM,RSECM,ILONM,ILMINM,RLSECM 
WRITE (6, 105) ILATM,IMINM,RSECM, ILONM, ILMINM , HLS ECM 
105 FORM AT (IX, 13, IX, 12, IX, F5. 2, IX, 14, IX, 12, IX, F5. 2) 

RMAST = ( (lABS (ILATM) * 60 + IMINH) ♦ 60 + RSECM) / 

1 RHOSEC 

IF (ILATM .LT. 0) RMAST = -RMAST 

RMASTL = ( (lABS (ILONM) * 60 + ILMINM) *60 +RLSECM) / 

1 RHOSEC 

IF (ILONM .GT. 0) RMASTL = TWOPI - RMASTL 
C 

C READ THE FIRST SLAVE STATION AND CHANGE THE LATITUDE 

C AND LONGITUDE INTO RADIANS. 
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c 

READ (5, 101) ILATS,IMINS,RSECS,rLONS,ILMINS,RLSECS, 

1 DELAY1 

WRITE (6,101) ILATS,IM INS , RSECS , I LONS , ILMIN S , RLSBCS , 

1 DELAY1 

101 FORaAT(1X,l3, 1X,I2,U ,F5. 2, 1X,I4, 1X,I2, 1X,F5.2,F9.2) 

RSLAV1 = ( (lABS (ILATS ) * 60 + IMINS) * 60 + RSECS) / 

1 RHOSEC 

IF (ILATS .LT. 0) RSLAV1 = -RSLAV1 

RSLAL1 = ( (lABS (ILONS) * 60 + ILMINS) =«= 60 + RLSECS) / 
1 RHOSEC 

IF (ILONS .GT. 0) RSLAL1 = TWOPI - RSLAL1 
C 

C READ THE SECOND SLAVE STATION AND CHANGE THE LATITUDE 

C AND LONGITUDE INTO RADIANS. 

C 

READ (5, 101) ILATS,IMI NS, RSECS, ILONS, ILMINS, RLSECS, 

1 DELAY2 

WRITE (6,101) IL ATS ,IM INS, RSECS, ILONS , ILMINS , RLSECS , 

1 DELAY2 

RSLAV2 = ( (lABS (ILATS) * 60 ^ IMINS) * 60 + RSECS) / 

1 RHOSEC 

IF (ILATS .LT. 0) RSLAV2 = -RSLA72 

RSLAL2 = ( (TABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 

IF (ILONS .GT. 0) RSLAL2 = TWOPI - RSLAL2 

WRITE (6,201) 

C 

C READ THE RECORDS POSITIONS AND OBSERVED RATES, THEN 

C COMPUTE THE TOTAL DISTANCE. 

C 

DO 550 I=1,IREC,1 

READ (5, 104) I LA T, I MIN , RSEC , I LON , I LMIN, RLSSC , 

1 HATE1,RATE2 



30 



104 FOEMAT(1X,I3,I3,F6.3, 15 ,13 , F6 . 3 , 3 X, 2F1 0 . 2) 

C 

C CONVERT LAT AND LONG TO RADIANS. 

C 

RPOST = ((lABS(ILAT) * 60 + IMIN) * 60 + RSEC) / 

1 RHOSEC 

IF(ILAT .LT. 0) RPOST = -RPOST 

RPOSTL = ((lABS(ILON) * 60 + ILMIN) * 60 + RLSEC) / 

1 RHOSEC 

IF(ILON .GT. 0) RPOSTL = TWOPI - RPOSTL 
C 

C COMPUTE DISTANCES AND AZIMUTHS FROM THE OBSERVED POINT 

C 

C MASTER 

C 

CALL INVER1 (A, F , RPOST , RPOSTL, RM AST, RMASTL, FASM , BAZM , 

1 DISTM) 

UDISTM = (RN * DISTM) / UNCOV 
CALL SECF AC (UDISTM, SFM) 

C 

C SLAVE1 

C 

CALL INVERT ( A , F , RPOST , RPO STL, RS L A V 1 , RSLAL1 , F AZ 1 , 3 AZ 1 , 

1 DIST1) 

UDIST1 = (RN * DIST1) / UNCOV 
CALL SECFAC (UDI ST1 ,SP 1) 

C 

C SLAVE2 

C 

CALL INVERT (A, F, RPOST , RPO ST L, RS L A V2 , RSL AL2 , F AZ 2 , 3AZ 2, 

1 DIST2) 

UDIST2 = (RN * DIST2) / UNCOV 
CALL SECFAC (UDIST2,SF 2) 

C 
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C COMPUTE THE RATES AND COMPARE TO THE OBSERVED RATES 

C 

TDM1 = UDIST1 - ODISTM + SF1 - SFM + DELAY 1 
DIFF1 = TDM1 - RAT El 
C 

TDM2 = UDIST2 - UDIST M + SF2 - SFM + DELAY2 
DIFF2 = TDH2 - RATE2 
C 

C COMPOTE THE LANE WIDTH IN METERS BASED ON EQUATION 

C 4.20 IN ELECTRONIC SURVEYING AND NAVIGATION - 

C LAURILA, PAGE 94. 

C 



BR1 = 


DABS (FAZ1 


- FAZM) 










3R2 = 


DABS (FAZ2 


- FAZM) 










WIDTHI 


= (DIFF1 


* UNCOV * 


0.5) 


/ 


DSIN (3R1 


* 0.5) 


WIDTH2 


= (DIFF2 


* UNCOV * 


0.5) 


/ 


DSIN (BR2 


* 0.5) 



C 

C WRITE THE POSITION OF VESSEL, COMPUTED RATE, OBSERVED 

C RATE, AND THE DIFF BETWEEN THEM. 

C 

WRITE (6,200) ILAT,IMI N,RSEC,ILON, ILMIN, RL3EC,RATS1, 
DIFF1 ,RATE2,DIFF2 

200 FORMAT( 1X,I3,I3 ,1X,F6 .3,15,13,1 X, F6. 3,F12. 2,F8 , 2, 

1 F12.2,F8.2) 

WRITE (6,201) 

C 

C XMEAN1 AND XMEAN2 ARE THE MEAN DIFFS BETWEEN THE 

C COMPUTED RATS AND THE OBERSERVED. XSLAV1 AND XSLAV2 

C ARE THE STORED DIFFS. 

C 

XMEAN1 = XMEAN1 + DIFF1 
XMEAN2 = XMEAN2 + DIFF2 
XMEAN3 = XMEAN3 + WIDTHI 
XMEAN4 = XMEAN4 + WIDTH2 
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550 

C 

C 

C 



C 



600 

C 



C 



C 



210 



211 



XSLAV1 (I) = DIFF1 

XSLAV2(I) = DIFF2 
CONTINUE 

COMPUTE THE MEAN AND STANDARD DEVIATION 

XMEAN1 = XMEAN1 / IREC 

XMEAN2 = XMEAN2 / IREC 

XMEAN3 = XMEAN3 / IREC 

XMEAN4 = XMEAN4 / IREC 

DO 600 I = 1 , IREC, 1 

VAR1 = VAR1 + ((XSLA7 1(I) - XMSAN1)**2) 

VAE2 = VAR2 + ((XSLA72(I) - XHEAN2)**2) 

VAR1 = VAR1 / (IREC - 1.0) 

VAR2 = VAR2 / (IREC - 1.0) 

SD1 = DSQRT(VARI) 

SD2 = DSQRT(VAR2) 

WRITE (6,201) 

WRITE (6,210) XMSAN 1,SD1 ,XMEAN3 
WRITS (6,201) 

WRITE(6,211) XMEAN2,S D2,XMEAN4 
WRITE (6,202) 

FORMAT(1X,' SLAVE #1, MEAN = »,F10.3, 

1 • STANDARD DEVIATION = •, 

2 F10.3,' DISTANCE IN METERS = ',F10.3) 
F0RMAT(1X,' SLAVE #2, MEAN = 

1 ' STANDARD DEVIATION = *, 

2 F10.3,* DISTANCE IN METERS = *,F10.3) 
STOP 

END 
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c 



SUBROUTINE SECFAC (UTD 1ST, SF) 

C 

C 

C THIS ROUTINE WILL COMPUTE THE SEA SECONDARY FACTOR 

C UTDIST = TOTAL DISTANCE 

C SF = SECONDARY FACTOR 

C 

IMPLICIT REAL*8 (A-H, 0-Z) 

C 

C 

C COEFFICIENTS 

C 

AO = 129.04323 

A1 = -0.40758 
A2 = 0.00064576813 

C 

BO = 2.741282 

B1 = -0.011402 
B2 = .00032774815 
C 

IF (UTDIST .GT. 537) GO TO 10 

SF = ( BO / UTDIST) + B1 + ( B2 * UTDIST) 

GO TO 20 

10 SF = ( AO / UTDIST) + A1 + ( A2 * UTDIST) 

C 

20 CONTINUE 
RETURN 
END 

C ========= = ====== = = ===^=== = ==:^= ====== = = === ====== = ^=== = ===^= 

SUBROUTINE INVERT (A, ? INV, GLAT 1 , GLON 1 , GL AT2 , GLO N2 , F AZ , 
1 BAZ,S) 

C 

C 
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C ♦♦♦ SOLUTION OF THE GEODETIC INVERSE PROBLSK AFTER 
C T.VINCSNTY MODIFIED RAINS FORD'S METHOD WITH HELMERT'S 

C ♦** ELLIPTICAL TERMS. EFFECTIVE IN ANY AZIMUTH AND AT 
C *** ANY DISTANCE SHORT OP ANTIPODAL STANDPOINT/FOREPOINT 
C MUST NOT BE THE GEOGRAPHIC POLE 

C 

C A IS THE SEMI-MAJOR AXIS OF THE REFERENCE ELLIPSOID 

C *** FINV IS THE FLATTENING (NOT RECIPROCAL) OF THE 
C *** REFERNECE ELLIPSOID LATITUDES AND LONGITUDES IN 
C *** RADIANS POSITIVE NORTH AND EAST FORWARD AZIMUTHS AT 
Q ««« both points returned in RADIANS FROM NORTH GEODESIC 
C *** DISTANCE S RETURNED IN UNITS OF SEMI-MAJOR AXIS A 
C 

C *** PROGRAMMED FOR CDC-66 00 BY LCDR L. PFEIFER NGS 
C **=<' ROCKVILLE MD 18FEB75. MODIFIED FOR IBM SYSTEM 360 
C by JOHN G GERGEH NGS ROCKVILLE MD 7507. 

C 

IMPLICIT REALMS (A-H,0-Z) 

DATA SPS/0.5D-1 3/,PI/3. 141 5926535898D0/ 

TWOPI=2.*PI 
R=1. -FINV 

TU1 = R*DSIN (GLAT 1) /DCOS (GLATl) 

' TU2=R*DSIN (GLAT2)/DCOS (GLAT2) 

CU1=1 ./DSQET (TU 1*TU1+ 1.) 

SU1=CU1#TU1 

CU2= 1 ./DSQRT (TU 2*TU2+ 1 . ) 

S = CU 1+CU2 
3AZ=S*TU2 
FAZ=BAZ*TU1 
X=GLON2-GLON1 
100 SX=DSIN(X) 

CX=DCOS (X) 

TU1=CU2*SX 

TU2=3AZ-SU1*CU2*CX 
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SY=DSQRT(TU1*TU 1+TU2*TU2) 

CY=S*CX+FAZ 
Y=DATAN2 (SY,CY) 

SA=S*SX/SY 

C2A=-SA*SA+1. 

CZ=FAZ+FAZ 

IF (C2A.GT.0.) CZ=-CZ/C2A+CY 
E=CZ*CZ*2.-1. 

C= ((-3.*C2A+4 .) +FINV+ 4.) *C2A*FI NV/1 6 . 

D = X 

X= ((E*CY*C+CZ) *SY*C+Y) ^S'SA 
X= (1 .-C) »X*FINV+GL0N2 -GL0N1 
IF (DABS (D-X) . GT.EPS) GOTO 100 
FAZ=DATAN2 (TD1,TU2) 

IF (FAZ.GE.TMOPI) FAZ= FAZ-T WOPI 
IF (FAZ.LT. 0. DO) FAZ= F AZ+T WOPI 
BAZ=DATAN2 (CU1*SX, BAZ *CX-SU1*C0 2) +PI 
IF (BAZ.GE.TWOPI) B AZ= BAZ-T WOPI 
IF (BAZ.LT.O.DO) BAZ=B AZ+T WOPI 
X = DSQRT ((1./R/R-1.) *C2A+1.) 

X= (X-2.)/X 
C=1. -X 

C= (X*X/4. +1 .) /C 
D= (0. 375»X*X- 1. ) *X 
X=E*CY 
S=1. -E-E 

S= ( ( ( (SY*SY*4.- 3.) *S*CZ*D/6.-X) *D/4. +CZ) *S Y*D + Y) *C*A*R 

RETURN 

END 

SENTRY 
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DATA SET EXAMPLE 



PEOGRAM LORAN 



SHIP'S POSITION OBSERVED LORAN RATES & ERRORS (E) 

LATITUDE LONGITUDE 9940-Y EY 9940-W EW 





(D-M 


-s) 




(D-M 


-S) 


(psec) (}isec) 


(;aS€C) (^S9C) 


36 


43 


45. 800 


121 


55 


27. 160 


42789.34 -0.49 


16294.04 -1.06 


36 


44 


3.400 


121 


55 


32. 340 


42791.13 -0.38 • 


16293.46 -1.10 


36 


44 


21. 180 


121 


55 


37. 390 


42793.04 -0.38 


16292. 73 -0. 99 


36 


44 


37. 490 


121 


55 


46. 950 


42795.13 -0.58 


16292.03 -1.06 


36 


44 


53. 260 


121 


55 


57.710 


42796.93 -0.51 


16291.43 -1.27 
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APPEN^X C 



PROG HAM LOPLC 



C PROGRAM LOPLC 

C 

C PROGRAM COMPOTES LINE OF POSITIONS FOR LORAN-C USING 

C THE SECONDARY FACTOR (SF) AND ADDITIONAL SECONDARY 

C FACTOR (ASF) THE SF AND ASF ARE BASED ON 

C SEMI-EMPIRICAL TD GRID. 

C 

IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION PHI (100) ,ELON (100) ,DISTM (100) ,DIST1 (100) , 

1 DIST2 (10 0) ,AXIS (13) , RF (13) ,AZ1 (100) , AZ2 (1 00) , AZM ( 1 00) 
DATA RHOSEC,PI, 0NCO7/2.062648062471D05, 

1 3. 1415926535898D0,299.693D0/ 

C 

CCC CCC**** *** ** **:^e****4e* 

DATA AXIS/6. 3782064D06, 6. 378388 D06 , 6. 377397 1 55D06 , 

1 6.37816 D06, 6. 37816D06, 6. 378249145006,6. 378165D06, 

2 6.378166006,6.378165 D06 , 6 . 378 1 4 5D0 6 , 6 . 3775634D06 , 

3 6.378245006,6.3781350006/ 

DATA RF/6. 3565838 D06, 2. 97 D 02, 2 . 99 15 28 12 85D02 , 

1 2.98 2500 2,2. 98 247 16 7 427D 02 , 2 . 9 3 46 5D02 , 2 .9 825D02 , 





2 2 


.983 D02,2. 983D02,2 .9825 


D02, 


6. 3562569D06 , 




3 2 


.983 D02,2. 98 26D02/ 






C** ♦♦♦ELLIPSOID OPTION NUMBER 






C 


1. 


CLARKE 1866 


8. 


MERCURY 


C 


2. 


INTERNATIONAL (HAY FORD) 


9. 


MARSHALL ISLAND 


C 


3. 


BESSEL 1841 


10. 


NAVY 8D 


C 


4. 


AND (AUSTRALIAN) 


11 . 


AIRY 


C 


5. 


1967 REFERENCE 


12. 


KRASSOWSKI 1940 


c 


6. 


CLARKE 1880 MOD 


13. 


WGS 1972 
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C 7. SAO 

C 

ccc ccc* ****** ******** ****** ************* 

C CC1-2 = ELLIPSOID NUMBER (K) 

C CC3-5 = NUMBER OF POINTS ALONG COAST (IREC) 

C 

READ (5, 100) K,IREC,ISTA1, ISTA2 

100 FORMAT(1X,I2,I4 ,213) 

TWOPI=2.*PI 

A = AXIS (K) 

F=1./RF (K) 

IF (F.LT.3.D-3) F=(A-1./F)/A 

C 

C READ MASTER STATION POSITIONS. 

C 

READ (5, 105) ILATH,IMI NM,RSECM, ILON M, ILMINM , RLSECM 
105 FORMAT ( 1X,I 3, 1X,I2,1X,F5.2,1X,I4,1X,I2,1X,F5.2) 

RMAST = ( (IA3S (ILATM) * 60 + IMINM) * 60 + RSECM) / 

1 R HO SEC 

IF (ILATM .LT. 0) RMAST = -RMAST 

RMASTL = ( (TABS (ILONM) * 60 + ILMINM) *60 +RLSECM) / 

1 RHOSEC 

IF (ILONM .GT. 0) RMASTL = TWOPI - RMASTL 
C 

C READ THE FIRST SLAVE STATION AND CHANGE THE LATITUDE 

C AND LONGITUDE INTO RA DI ANS . DELAY IS THE CODING DELAY 

C AND BIAS IS THE OFFSET IN MICROSECONDS. 

C 

READ (5, 101) ILATS,IMI NS ,R SECS , I LON S, ILMINS , RLSECS , 

1 DELAY1,3IAS1 

101 FORM AT ( 1X,I3, 1X,I2,1X,F5.2,1X,I4,1X,I2,1X,F5.2, 

1 F9.2,F7.3) 

RSLAV1 = ( (IA3S (ILATS ) * 60 + IMINS) * 60 RSECS) / 
1 RHOSEC 
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-RSLAV1 



IF(ILATS .LT. 0) RSLA71 = 

RS1AL1 = ( (lABS (ILONS ) ♦ 60 + ILMINS) * 60 + RLSECS) / 
1 RKOSEC 

IF (ILONS .GT. 0) RSLAL1 = TWOPI - RSLAL1 
C 

C READ THE SECOND SLAVE STATION AND CHANGE THE LATITUDE 

C AND LONGITUDE INTO RADIANS. 

C 

READ (5, 101) HATS, IMI NS, RSECS, ILONS, ILHINS, RLSECS, 

1 DELAY2,BIAS2 

RSLAV2 = ( (TABS (ILATS) * 60 + IMINS) * 60 + RSECS) / 

1 RHOSEC 

IF (ILATS .LT, 0) RSLAV2 = -RSLAV2 

RSLAL2 = ( (TABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 

IF (ILONS .GT. 0) RSLAL2 = TWOPI - RSLAL2 
C 

C READ COAST POINT LAT AND LONG AND CONVERT TO RADIANS. 

C STORE LAT IN PHI AND LONG IN ELON. 

C 

DO 500 I=1,IR2C,1 

READ (5,102) INUM,ILAT,IMIN,RSEC,ILON,ILMIN, RLSEC 
102 FORM AT ( 1 X, 14, 14, 13, P7 .3, 1 5,13, F 7. 3) 

PHI (I) = ((lABS(ILAT) * 60 + IMIN) * 60 + RSEC) / 

1 RHOSEC 

IF(ILAT .LT. 0) PHI (I) = -PHI (I) 

ELON (I) = ((lABS(ILON) * 60 + ILMIN) * 60 + RLSEC) / 

1 RHOSEC 

IF(ILON .GT. 0) ELON (I) = TWOPI - ELON (I) 

C 

C COMPUTE DISTANCE FROM COAST POINT TO MASTER AND 

C SLAVE STATIONS. SET UP COMPUTATION AND CALL INVERl 

C 

PI = PHI (I) 
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El = ELON(I) 

CALL INVER1 (A, F, PI, El ,RMAST , RMASTL , AZF , AZB , S) 
DISTM(I) = S 
AZM(I) = AZF 

CALL INVER 1 (A,F,P1 ,E1 , RSL A V 1 , RS LAL 1 , AZF,AZB,S) 

DIST1 (I) = S 
AZ1(I) = AZF 

CALL INVERT (A, F, PI, El ,RSLAV2, RSLAL2 , AZF , AZB , S) 

DIST2 (I) = S 
AZ2(I) = AZF 

WRITE (6,210) I, INUM,I LAT,i:iIN, RSEC,ILCN,ILi3IN, RLSEC, 
1 AZM (I) ,AZ1 (I) , AZ2 (I) 

210 FORM AT( 1X,I4,I4 ,14 ,13 , F7. 3 , 15 , 1 3 , F7 . 3 , 3F 1 0 . 3 ) 

500 CONTINUE 
C 

C READ THE NUMBER OF RECORDS AND TKEIR POSITIONS. 

C COMPUTE THE TOTAL DISTANCE OVER LAND AND SEA. 

C TDISTM = TOTAL DISTANCE FROM MASTER TO DATA POINT. 

C TDISTS = TOTAL DISTANCE FROM SLAVE TO DATA POINT. 



C 



C 



COMPUTE THE FORWARD AND BACK AZIMUTHS 



C 



FAZM AND FAZS = FORWARD AZIMUTH TO THE MASTER AND 



C 



SLAVE STATION 



C 



EAZM AND 3AZS = BACK AZIMUTH TO THE MASTER AND 



C 



SLAVE STATION 



C 



READ (5, 103) JREC 
103 FORM AT (IX, 14) 

C 



DO 550 1=1, JREC, 1 



READ (5, 104) ILAT ,IMIN, HSEC ,ILON, ILMIN , RLSEC, HATH 1 , RATE2 
104 FORM AT (IX, 13,13 ,F6. 3, 15,1 3, F6. 3, F 9.3, F9. 3) 



C 

C CONVERT LAT AND LONG TO RADIANS. 
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c 

RPOST = ((lABS(ILAT) ♦ 60 + IMIN) * 60 + RSEC) / 

1 RHOSEC 

IF(ILAT .LT. 0) RPOSI = -RPOST 

RPOSTL = {(lABS(ILON) * 60 + ILMIN) * 60 + RLSEC) / 

1 RHOSEC 

IF(ILON .GT. 0) RPOSTL = TWOPI - RPOSTL 
C 

C COMPUTE DISTANCES AND AZIMUTHS FROM THE OBSERVED 

C POINT TO THE INTERPOLATED SHORE POINT AND WRITE. 

C 

ISTA = 1 

CALL SHORP? (IREC,RPOS T, RPOSTL , R MAS T , RMASTL , UTDI SM, 

1 USDIST,ULDIST, FAZM, B AZM) 

CALL SECFAC (DTD ISM , US DIST , ULDI5 T , B AZM , I STA , SFM ) 

C 

CALL SHORPT (I REC, RPOS T, RPOSTL , RSLA VI , RSLAL 1 , UTDIS 1 , 

1 USDIST,ULDIST,FAZ1,BAZ1) 

CALL SECFAC (UTDIS 1 , OS DIST , OLDIS T , BAZ 1 , ISTA 1 , SF 1 ) 

C 

CALL SHORPT (IREC,R?OS T, RPOSTL, RSLAV2, RSLAL2 , UTDIS2, 

1 USDIST,OLDIST, FAZ2,3 AZ2) 

CALL SECFAC (UTD IS2 ,US DIST , ULDIS T , BAZ2 , ISTA2 , SF2) 

C 

C COMPUTE THE HATES AND COMPARE TO THE OBSERVED RATES 

C 

TDM1 = OTDIS1 - UTDISM + SF1 - SFM + DELAYl + BIAS1 
DIFF1 = TDMI - RATS1 
C 

TDM2 = UTDIS2 - UTDISM + SF2 - SFM + DELAY2 *■ BIAS2 
DIFF2 = TDM2 - RAT S2 
C 

C COMPUTE THE LANE WIDTH IN METERS BASED ON EQUATION 

C 4.20 IN ELECTRONIC SURVEYING AND NAVIGATION - LAURILA, 
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PAGE 94. 








BRI = DABS (FAZ1 


- FAZM) 






IF(BR1 .GT. PI) 


BRI = TWO PI - 


BRI 




BR2 = DABS (FAZ2 


- FAZ M) 






IF(BR2 .GT. PI) 


BR2 = TWOPI - 


3R2 




WIDTH1 = (DIFF1 


♦ UNCOV * 0.5) 


/ DSIN(BR1 ♦ 


0. 5) 


WIDTH2 = (DIFF2 


* UNCOV * 0.5) 


/ DSIN(BR2 * 


0. 5) 



c 

C WRITS THE POSITION OF VESSEL, OBSERVED RATES, 

C AND THE DIFFERENCES BETWEEN THEM. 

C 

WRITE (6 ,200) ILAT,IMI N, RSEC, ILON , ILMIN , RLSSC , R ATE 1 , 

1 DIFF1, RATE2, DIFF2 

200 FORM AT (IX, 1 2, 13 ,F7 . 3 , 14 ,I 3 , F7 . 3 , F 1 0 . 2 , F6 . 2 , F 10 . 2 , F6 . 2) 

C 

C XMEAN1 AND XHEAN2 ARE THE MEAN DIFFS BETWEEN THE 

C COMPUTED RATE AND THE OBERSSRVBD. XSLAV1 AND XSLAV2 

C ARE THE STORED DIFFS. 

C 

XMSAN1 = XMEAN1 + DIF FI 

XMEAN2 = XMEAN2 + DIFF2 

XMEAN3 = XMSAN3 + WIDTH 1 

XMEAN4 = XMEAN4 + WIDTH2 

XSLAVI(I) = DIFF1 
XSLAV2(I) = DIFF2 
550 CONTINUE 
C 

C COMPUTE THE MEAN AND STANDARD DEVIATION 

C 

XMEAN1 = XMEAN1 / JREC 

XMSAN2 = XMEAN2 / JREC 

XMEAN3 = XMEAN3 / JREC 

XMEAN4 = XMEAN4 / JREC 
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c 

DO 600 I = 1, JSEC, 1 

VAR1 = VAB1 + ((XSLA7 1(I) - XMEAN1)**2) 

600 VAR2 = 7AR2 + ((XSLA72(I) - XHEAN2) **2) 

C 

7AR1 = 7AR1 / (JREC - 1.0) 

7AR2 = 7AR2 / (JREC - 1.0) 

C 

SD1 = DSQRT (7AR 1) 

SD2 = DSQRT (7AR2) 

C 

WRITE (6,201) 

WRITE (6,210) XMEAN 1,S D1 ,XMEAN3 
WRITE (6,201) 

WRITE (6,211) XHEAN2,S D2,XMEANa 
WRITE (6,202) 

210 ?ORMAT(1X,' SLA7E #1, MEAN = ’,F10.3, 

1 ' STANDARD DE7IATI0N = ', 

2 F10.3,’ DISTANCE IN METERS = «,?10.3) 

211 FORM AT (IX,’ SLA7E #2, MEAN = ',F10.3, 

1 ' STANDARD DS7IATI0N = ', 

2 F10.3,’ DISTANCE IN METERS = ',?10.3) 

WRITE (6,201) 

201 FORMAT (ia ) 

202 FORMAT(IHI) 

STOP 

END 

C== ================================== = =========== ============ 

S03RO0TINE SHORPT (IRE C, RPOST , RPOSTL, RCOIIT, RCONTL , ODIST , 
1 OSD, OLD, FA, 3A) 

C 

C 

C SOBROOTINE SHORPT WILL SELECT A POINT ALONG THE SHORE 

C WHICH IS OUTLINED FROM NORTH BY SELECTED POINTS WITH 
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C KNOWN LATITUDES AND LONGITUDES. THE SHORE POINT IS 

C INTERPOLATED BETWEEN TWO KNOWN POINTS USING THE TOTAL 

C DISTANCE BETWEEN THE POSITION AND THE CONTROL STATION 

C AND THE AZIMUTH BETWEEN THE SHORE POINTS AND THE 

C RECEIVERS POSIT lON- 

C 

IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION FAZM(IOO) 

COMMON/SHORE/PHI (100) ,ELON (100) , UNCOV,RN, A, F 
DATA PI/3, 14159 26535398D0/ 

C 

TWOPI = 2.0 * PI 
C 

CALL INVER1 (A,F,RPOST , RPO STL, RCONT, RCONTL, FA , B A, RDISTT) 
DDIST = (RN * RDISTT) / UNCOV 
C 

RCOMP = 99999.99 
C 

DO 10 J=1 ,IREC, 1 
PI = PHI(J) 

El = ELON(J) 

CALL INVERT (A, F, PI, El , RCONT , 2C0 NTL , AZF , AZ3 , R DI STL) 

CALL INVERT (A ,F , RPOST , RPOSTL, PI , E 1 , FAZ , BAZ , RDISTS) 

OSD = (RN * RDISTS) / UNCOV 
OLD = (RN * HDISTL) / DNCOV 
FAZM (J) = AZE 
C 

C COMPUTE THE DIFFERENCE BETWEEN THE TOTAL DISTANCE 

C (DDIST) AND THE SUMMATION OF THE DISTANCE OVER THE WATER 

C (USD) AND THE DISTANCE OVER THE LAND (OLD) . IT THE 

C DISTANCE IS LESS THAN RCOMP, UPDAT RCOMP AND JSTA. 

C JSTA IS THE CLOSEST POINT ALONG THE SHORELINE WHICH IS 

C NEAR THE EM PROPAGATION PATH. 

C 
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FDIF? = DABS(UDIST - (USD + OLD)) 

IF(FDIFF .GT. RCOMP) GO TO 10 
RCOMP = FDIFF 
JSTA = J 
10 CONTINUE 

C 

C NOW DETERMINE THE INTERPOLATED LATITUDE AND LONGITUDE 

C SHORE POINT USING AZIMUTH PERCENTAGE. 

C 

lUPPER = JSTA - 1 
ILOWER = JSTA + 1 
AZMU = FAZM(IUPPER) 

AZML = FAZa (ILOWER) 

IF (AZMU .GT. FA .AND. FA .GE. FAZM (JSTA) ) ICH = lUPPEH 
IF (FAZM (JSTA) . GE. FA .AND. FA .GT. AZML) ICH = ILOWER 
C 

RADJ = 1.00 - DABS ((FAZM(ICH) - FA) / (FAZM (ICH) - 
1 FAZM (JSTA) ) ) 

C 

RNWPHI = PHI (JSTA) + ((PHI (ICH) - PHI (JSTA) ) * RADJ) 

IF (ELON (JSTA) . LE. ELON(ICH)) RADJ = 1.0 - RADJ 
RNWELN = ELON (JSTA) + ((ELON (ICH) - ELON (JSTA)) * RADJ) 

C 

CALL TODHS (RNWPHI, IDG , MIN, SEC) 

RHOLD = TWOPI - RNWELN 

CALL TODMS (RHOL D,IDG L , MIN L, SECL) 

C 

CALL INVERT (A, F ,RNWPH I, RNWELN, RCONT, RCONTL , AZ, BZ , RDISTL ) 

CALL INVERT ( A , F , RPOST , RPO STL , RN WP HI , RNWELN , A Z, 5Z , RDISTS ) 

USD = (RN * RDISTS) / UNCOV 

ULD = (RN * RDISTL) / UNCOV 

RETURN 

END 

C== ================== ========== ============================ 
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SUBROUTINE SECF AC (UTD 1ST, USDIST , ULDIST, AZI , I STA , SF) 



C- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



THIS ROUTINE WILL COMPUTE THE LAND/SEA SECONDARY FACTOR 
UTDIST = TOTAL DISTANCE 

USDIST = DISTANCE OVER THE SEA WATER PATH 
OLDIST = DISTANCE OVER THE LAND PATH 
AZI = AZIMUTH FROM NORTH. 



MASTER 

W 

X 

Y 



1 (ISTA) 

2 
3 



IMPLICIT REAL*8 (A-H,0-Z) 





SI 


= (. 795 / 


USDIST) + 0. 439 + 




(.00245 * 


USDIST) 


c 


















IF 


(UTDIST .GT 


. 540) GO TO 10 












S2 


= (3.188 / 


UTDIST) - 0.594 


+ 


(.000329 




UTDIST) 




GO 


TO 20 












10 


S2 


= (128.8 / 


UTDIST) + 0. 187 




(.000652 


•X, 


UTDIST) 


c 
















20 


IF 


(USDIST .GT 


. 540) GO TO 30 












S3 


= (3.188 / 


USDIST) - 0,594 




(.000329 




USDIST) 




GO 


TO 40 












30 


S3 


= (128.8 / 


USDIST) + 0. 187 




(.000652 




USDIST) 


c 
















40 


IF 


(ULDIST .GT 


. 540) GO TO 50 












S4 


= (3.188 / 


ULDIST) - 0.594 




(.000329 


4: 


ULDIST) 




GO 


TO 60 












50 


S4 


= (128.8/ 


ULDIST) + 0.187 




(.000652 




ULDIST) 


C 
















60 


S5 


= 1.42 8 + 


(. 001 58 * UTDIST) 












S6 


= 1. 428 + 


(. 00153 * ULDIST) 
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70 



80 

90 

100 

C 

200 



500 



C 



C 



TAZI = 2. 0 ♦ AZI 
GOTO (70,80, 90 ,100) , ISTA 
RHOLD = (1. 01 0*DSIN (AZI) ) 

1 - (.893*DSIN (TAZI) ) 



GO TO 200 

RHOLD = (. 323«DCOS ( AZI) ) 

GO TO 200 

RHOLD = (.942*DCOS (TAZI) ) 

GO TO 200 

RHOLD = (.588*DSIN (TAZI) ) 



(. 196 + DCOS (AZI) ) 
(. 355*DCOS (TAZI) ) 

(.711 *DSIN (TAZI) ) 



55 = S5 + RHOLD 

56 = S6 + RHOLD 

SF = 0.5 * (S5 + S6 - Si + S2 + S3 - SU) 
WRITE (6,5 00) SI ,S2,S3 ,S4,S4,S5, S6 ,S? 
FORMAT (IX, 7F15. 5) 

RETURN 

END 



SUBROUTINE INVER1 ( A, F INV, GLAT1 , GLON 1 ,GLAT2 , GLON2 ,FAZ, 
1 3AZ,S) 



See Appendix B for subroutine INVERI. 
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Short points around Monterey Bay, California. The point 
are used to interpolate geodetic point for computarion 
of the distance over land and the distance over sea. 



NO. LATITDDB LONGITUDE 



1 


36 


57 


18 . 


.606 


122 


05 


37 , 


,525 


2 


36 


56 


59 . 


.264 


122 


03 


01 , 


,817 


3 


36 


57 


05 . 


.076 


122 


01 


31 , 


,701 


4 


36 


57 


49 . 


.538 


122 


01 


07 , 


857 


5 


36 


57 


17 . 


.949 


121 


58 


19 . 


830 


6 


36 


58 


08 . 


.589 


121 


57 


07 , 


,288 


7 


36 


58 


32 . 


.140 


121 


55 


10 . 


, 083 


8 


36 


58 


01 , 


.498 


121 


53 


57 , 


,390 


9 


36 


56 


46 . 


.115 


121 


52 


22 , 


313 


10 


36 


55 


38 . 


. 140 


121 


51 


24 , 


399 


11 


36 


53 


13 . 


.806 


121 


49 


46 , 


,743 


12 


36 


49 


38 . 


.384 


121 


47 


48 , 


,8 95 


13 


36 


47 


39 . 


.241 


121 


47 


10 , 


,818 


14 


36 


46 


27 . 


.554 


121 


47 


39 . 


,637 


15 


36 


44 


56 . 


.717 


121 


47 


52 , 


,4 16 


16 


36 


41 


14 . 


.439 


121 


48 


32 , 


, 642 


17 


36 


39 


17 . 


.211 


121 


49 


28 , 


,533 


18 


36 


37 


31 . 


.128 


121 


50 


31 , 


,723 


19 


36 


36 


23 . 


.446 


121 


51 


34 , 


, 833 


20 


36 


36 


03 . 


.628 


121 


52 


5 0 , 


879 


21 


36 


36 


24 , 


.782 


121 


53 


48 , 


,453 


22 


36 


37 


18 . 


.151 


121 


54 


11 . 


628 


23 


36 


38 


00 . 


.300 


121 


55 


57 . 


,538 
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DATA SET EXAMPLE 



PROGRAM LOPLC 



SHIP'S POSITION OBSERVED LORAN RATES & ERRORS (E) 

LATITUDE LONGITUDE 9940-Y EY 9940-W EW 





(D-M 


-S) 




(D-M 


-S) 


ijisec) 


(fisec) 


^sec) 


(^sec) 


36 


43 


45. 800 


121 


55 


27. 160 


42789. 34 


1 . 14 


16294.04 


1. 32 


36 


44 


3.400 


121 


55 


32. 340 


42791.13 


1 . 25 


16293 .46 


1.29 


36 


44 


21. 180 


121 


55 


37. 390 


42793.04 


1.25 


16292.73 


1. 40 


36 


44 


37.490 


121 


55 


46.950 


42795. 13 


1.05 


16292.03 


1.33 


36 


44 


53. 260 


121 


55 


57. 710 


42796.93 


1 . 12 


16291 .43 


1. 11 
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APP ENDI X D 



PROGRAM LORTAB 



C PROGRAM LORTAB 

C 

C PROGRAM COMPOTES LINE OF POSITION FOR LORAN-C USING 

C SF SALT WATER CORRECTION FACTOR AND DMAHTC CALCULATED 

C OR FIELD OBSERVED ASF CORRECTIONS. 

C 

C 

IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION AXIS ( 13) ,RF (13) 

DIMENSION XSLAV1 (1000) ,XSLAV2 (1000) 

DATA RHOSEC,PI, UNCOV, RN/2 . 06 264 80 624 71 DOS , 

1 3.1415926535898D0,299.792458D0 ,1 .00033800/ 

DATA XMEAN1,XMEAN2,VAR1//AH2/0. OODO, O.OODO, 

1 O.OODO, 00. ODO/ 

DATA XMEAN3,XMEAN4/0. OODO, O.OODO/ 

C 

CCC CC C** ******** 

DATA AXIS/6. 3782064D0 6,6. 373388 D06 , 6. 3773971 55D06 , 

1 6.37816 D06, 6. 37816006,6. 378249145006,6. 373165D06, 

2 6.378166006,6.378165 D06 , 6. 378 1 45D06, 6. 3775634D06, 

3 6.378245006,6.3781350006/ 

DATA RF/6. 3565838 D06, 2.97 D02 , 2 .99 1528 12 85D02 , 

1 2.98 2500 2,2. 98 247 16 7 4 27D 02, 2.93 465002,2.9 8 25 DO 2, 

2 2.983 D02,2. 983D02,2 .9825 D02, 6. 3562569D06 , 

3 2.983 002,2.9826002/ 

C*****SLLI?SOID OPTION NUMBER 
C 1. CLARKE 1866 

C 2. INTERNATIONAL (HA I FORD) 

3. BESSEL 1841 



C 
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8- MERCURY 
9. MARSHALL ISLAND 
10. NAVY 8D 



c 


4. AND (AUSTRALIAN) 


11 . 


AIRY 


c 


5. 1967 REFERENCE 


12. 


KRASSOWSKI 1940 


c 


6. CLARKE 1880 MOD 


13. 


WGS 1972 


c 


7. SAO 







c 

CCC CCC* ****** ♦♦ ***♦ ♦* * 

C CCl-2 = ELLIPSOID NUMBER (K) 

C CC3-5 = NUMBER OF POINTS ALONG COAST (IREC) 

C 

READ (5, 100) K,IREC,ID 1,ID2 
100 FORMAT (IX, 12, 14 ,213) 

TWOPI=2.*PI 

A=AXIS(K) 

F=1./RF (K) 

IF (F.LT.3. D-3) F=(A-1./F)/A 

C 

C READ MASTER AND SLAVE STATIONS POSITIONS 

C THE FIRST RECORD IS THE NUMBER OF MASTER AND SLAVS 

C STATIONS FOR THE PARTICULAR CHAIN. 

C 

WRITE (6,202) 

202 FORMAT(IHI) 

WRITE (6,201). 

201 F0RMAT(1H ) 

READ (5, 105) ILATM, IMI NM,RSECM, ILONM, ILMINM , RLSECM 
WRITE (6 ,105) ILATM ,IMINM,RSECM, ILONM, ILMINM, RLSECM 
105 FORMAT(1X,I3, 1X,I2,1X ,F5. 2, 1X,I4, 1X,I2, 1X,F5.2) 

RMAST = ( (lABS (ILATM) * 60 + IMINM) * 60 + RSECM) / 

1 RHOSEC 

IF (ILATM .LT. 0) RMAST = -RMAST 

RMASTL = ( (lAES (ILONM) * 60 + ILMINM) *60 +RLSECM) / 
1 RHOSEC 

IF (ILONM .GT. 0) RMASTL = TWOPI - RMASTL 
C 
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C READ THE FIRST SLAVE STATION AND CHANGE THE LATITUDE 

C AND LONGITUDE INTO RADIANS. 

C 

READ (5, 101) ILATS, IMI NS ,RSECS , I LO NS , ILMINS , RLS ECS , 

1 DELAY1 

WRITE (6,101) ILATS ,IM INS,RSECS, ILONS , ILMINS, 5LSECS, 

1 DELAY1 

101 FORMAT (IX, 13, IX , 12 , 1 X , F5 . 2 , 1 X , I 4 , 1 X, 1 2 , IX , F5 . 2 , F9 . 2) 

RSLA71 = ( (lABS (ILATS) ♦ 60 + IHINS) * 60 + RSECS) / 

1 RHOSEC 

IF (ILATS .LT. 0) RSLAV1 = -RSLA71 

RSLAL1 = ( (lABS (ILONS) * 60 + ILMINS) ♦ 60 + RLSECS) / 
1 RHOSEC 

IF (ILONS .GT. 0) RSLAL1 = TWOPI - R3LAL1 
C 

C READ THE SECOND SLAVE STATION AND CHANGE THE LATITUDE 

C AND LONGITUDE INTO RADIANS. 

C . 

READ (5, 101) ILATS,IMINS, RSECS, ILONS, ILMINS, RLSECS, 

1 DELAY2 

WRITE (6,101) IL ATS, IM INS, RSECS, ILONS, ILMINS, RLSECS, 

1 DELAY2 

RSLA72 = ( (lABS (ILATS ) * 60 + IHINS) * 60 + RSECS) / 

1 RHOSEC 

IF (ILATS .LT. 0) RSLAV2 = -RSLA72 

RSLAL2 = ( (TABS (ILONS) ♦ 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 

IF (ILONS .GT. 0) RSLAL2 = TWOPI - RSLAL2 

WRITE (6,201) 

C 

C READ THE RECORDS POSITIONS AND OBSERVED RATES, THEN 

C COMPUTE THE TOTAL DISTANCE. 

C 

DO 550 I=1,IR3C,1 



103 



READ (5, 104) ILAT,IMIN , RSSC , ILON , ILMIN , RLSEC , 

1 RATE1,RATE2 

104 FORMAT (IX, 13,13 ,F6. 3, 15 ,I 3 , F6 . 3 , 3X, 2F 1 0 . 2) 

C 

C CONVERT LAT AND LONG TO RADIANS. 

C 

RPOST = ((lABS(ILAT) * 60 + IHIN) * 60 + RSEC) / 

1 RHOSSC 

IF (ILAT .LT. 0) RPOST = -RPOST 

HPOSTL = ((IA3S(IL0N) * 60 + ILMIN) * 60 + RLSEC) / 

1 RHOSEC 

IF(ILON .GT. 0) RPOSTL = TWOPI - RPOSTL 
C 

C COMPUTE DISTANCES AND AZIMUTHS FROM THE OBSERVED POINT 

C 

C MASTER 

C 

CALL INVERT (A ,F , RPOST , RPOSTL, RM AS T , RMASTL, F AZM , 3 A2M , 

1 DISTM) 

UDISTM = (RN ♦ DISTM) / UNCOV 
CALL SECFAC (DDISTM,SF M) 

C 

C SLAVS1 

c 

CALL INVERT (A,f , RPOST , RPO STL, RS LA VI , RSLAL 1 , FAZI , BAZ1, 

1 DIST1) 

UDIST1 = (RN ♦ DIST1) / UNCOV 
CALL SECFAC (UDIST1 ,SF 1) 

C 

C SLAVE2 

C 

CALL INVERT (A , F , RPOST , RPO STL, RS L A V2 , RSL AL2 , F AZ2 , 3 AZ 2 , 

1 DIST2) 

UDIST2 = (RN ♦ DIST2) / UNCOV 
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CALL SSCFAC{aDIST2,SF2) 

C 

C DETERMINE THE ADDITIONAL SECONDARY CORRECTORS FROM THE 

C LORAN-C CORRECTION TABLE FOR THE WEST COAST CHAIN 9940 

C 

CALL TABLE(RPOST,R?OSTL,ID1,ASF1) 

CALL TABLE(RPOST,RPOSTL,ID2,ASF2) 

C 

C COMPUTE THE RATES AND COMPARE TO THE OBSERVED RATES 

C 

TDM1 = UDIST1 - UDISTM + SF1 - SFM + ASF1 + DELAY1 
DIFF1 = TDM1 - RATE1 
C 

TDM2 = UDIST2 - UDISTM + SF2 - SFM + ASF2 + DELAY2 
DIFF2 = TDM2 - RATE2 
C 
C 

C COMPUTE THE LANE WIDTH IN METERS BASED ON EQUATION 

C 4.20 IN ELECTRONIC SURVEYING AND NAVIGATION - 

C LAURILA, PAGE 94. 

C 



BE1 = DABS (FAZ1 


- FAZM) 










3R2 = DABS(FAZ2 


- FAZ M) 










WIDTH1 = (DIFF1 


* UNCOV * 


0.5) 


/ 


DSIN(BR1 * 


0. 5) 


WIDTH2 = (DIFF2 


♦ UNCOV * 


0.5) 


/ 


DSIN(BR2 * 


0.5) 



C 

C WRITE THE POSITION OF VESSEL, COMPUTED RATE, OBSERVED 

C RATE, AND THE DIFF BETWEEN THEM. 

C 

WRITE(6,200) ILAT,IMI N,RSEC,IL0N,ILMIN, RLSEC,RAT31, 
DIFF 1 ,RATE2, DIF F2 

200 F0RMAT(1X,I3,I3,1X,F6 . 3 ,1 5 , 13 , 1 X , F6. 3 , F 12. 2 , F8 . 2 , 

1 F12.2,F8.2) 

WRITE (6,201) 
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c 

c 

c 

c 

c 



550 

c 

c 

c 



c 



600 

c 



c 



c 



XMEAN1 AND XMEAN2 ARE THE MEAN DIFFS BETWEEN THE 
COMPUTED RATE AND THE OBERSERVED. XSLAV1 AND XSLAV2 
ARE THE STORED DIFFS. 

XMEAN1 = XMEAN1 + DIFF1 
XMEAN2 = XMEAN2 + DIFF2 
XMEAN3 = XMEAN3 + WIDTH1 
XMEAN4 = XMEAN4 + WIDTH2 
XSLAV1 (I) = DIFF1 

XSLA?2(I) = DIFF2 
CONTINUE 

COMPUTE THE MEAN AND STANDARD DEVIATION 

XMEAN1 = XMEAN1 / IREC 
XMEAN2 =-XMEAN2 / IREC 
XKEAN3 = XMEAN3 / IREC 
XMEAN4 = XMEAN4 / IREC 

DO 600 I = 1 , IREC, 1 

VAR1 = VAR1 + ((XSLA7 1(I) - XMEAN1)^‘*2) 

VAR2 = 7AR2 + ((XSLA72(I) - XMEAN2)**2) 

VAR1 = VAR1 / (IREC - 1.0) 

VAR2 = VAR2 / (IREC - 1.0) 

SD1 = DSQRT (VAR 1) 

SD2 = DSQRT (7AR2) 

WRITE (6,201) 

WRITE (6,210) XMEAN1,S D1 ,XMEAN3 
WRITE (6 ,201) 

WRITE (6,211) XMEAN2,S D2,XMEAN4 
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210 



• ,F10.3, 



WRITS (6,202) 

FORMAT(1X,' SLAVE #1, MEAN = 

1 • STANDARD DEVIATION = 

2 F10.3,» DISTANCE IN METERS = »,F10.3) 

211 FORMAT(1X,» SLAVE #2, MEAN = »,F10.3, 

1 • STANDARD DEVIATION = 

2 F10,3,» DISTANCE IN METERS = ',?10.3) 

STOP 

C================ = === ====== ============ ========= = === = = == 

SUBROUTINE SECF AC ( UTD 1ST, SF) 

C 

C 

C THIS ROUTINE WILL COMPUTE THE SEA SECONDARY FACTOR 

C UTDIST = TOTAL DISTANCE 

C SF = SECONDARY FACTOR 

C 

IMPLICIT REAL*8 (A-H,0-Z) 

C 

C 

C COEFFICIENTS 

C 



AO 


= 


129.04323 












A1 


= 


-0.40758 












A 2 


= 


0.0006457 


6813 










BO 




2.741282 












B1 




-0.011402 












B2 




.000327748 


15 










I? 


(UTDIST .GT. 


537) 


GO TO 


10 






SF 




( BO / DTD 


1ST) 


+ B1 + 


( 


32 


* UTDIST) 


GO 


TO 


20 












SF 




( AO / UTD 


1ST) 


+ A1 + 


( 


A2 


* UTDIST) 



C 



107 



20 



CONTINUE 

BETURN 

END 



C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 



c 



SUBROUTINE TASLE(RLAT ,RLON,ID, ASF) 



SUBROUTINE TABLE SELECTS THE PBOPER ASF CORRECTOR FROM 
THE LORAN-C CORRECTION TABLE PUBLISHED BY THE DEFENSE 
MAPPING AGENCY. 

RLAT = POSITION LATITUDE IN SECONDS 
RLON = POSITION LONGITUDE IN SECONDS 
ID = LORAN-C CHAIN IDENTIFIER 
W = 1 
X = 2 
Y = 3 

ASF = ADDITIONAL SECONDARY FACTORS 



THE FOLLOWING TABLES OF ASF CORRECTORS ARE FOR 
MONTEREY BAY, CALIFORNIA - 9940 -H, -X, -Y. 



IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION TA3LEW(3,5) ,TAB LEX (3 , 5) , T ABLEY (3 , 5) 



DATA TABLEW/ 

1 

2 

3 

4 



1. 6D0, 1, 6D0, O.ODO, 
1. 6D0, 1- 5D0, 1 . 4D0, 
1. 5D0,1. 4D0,1 .6D0, 
1. 4D0,1. 3D0,1 .5D0, 
1. 3D0,0- ODO, O.ODO/ 



DATA TABLSX/ -0 . 9D0, - 0. 9D0 , O.ODO, 
-0 .9D0,- 0.9D0,-0.8D0 , 
-1 .ODO,- 1. 0D0,-1 .ODO , 
-1 .ODO,- 1.2D0,-1.1D0, 



1 

2 

3 
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4 -1.1D0, O.ODO, O.ODO/ 

C 

DATA TABLE!/ 0. 2D0 , 0. 2D0, 0 . ODO, 

1 0. 3D0,0. 3D0,0, 3D0, 

2 0. 3D0,0. 2D0,0. 4D0, 

3 0, 4D0 ,0. 3D0,0.6D0, 

4 0. 5D0,0. ODO, O.ODO/ 

DATA RHOS EC, PI/ 2. 06264306 247 1D05, 3. 14 159265358 98 DO/ 

C 

TWOPI = PI * 2.0 
C 

C CONVER RLAT AND RLON TO SECONDS 

C 

HLAT = RLAT ♦ RHOS EC 
HLON = TWOPI - RLON 
HLON = HLON * RHOS EC 
C 

C STARTING LAT AND LONG FOR SEARCH 

C LAT = 37/00/00.0 LONG = 122/05/00.0 

C 

c 

C DETERMINE THE ASF CORRECTOR FOR THE LORAN-C COMBINATION 

C 

C LATITUDE 

C 

SLAT = 133200.0 
SLON = 439500-0 
RMID = 300.0 
RDIFF = 150-0 
C 

DO 10 J=1, 5, 1 
SLAT = SLAT - RMID 
ULAT = SLAT + RDIFF 
VLAT = SLAT - RDIFF 
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i 



% 




TO 


IF (hLAT .LT. ULAT .AND. 


hLAT .GE. 


VLAT) 


GO 


TO 


T5 


T5 


CONTINUE 












C 














C 


LONGITUDE 












c 


DO 30 I=T,3,T 














SLON = SLON - RaiD 














ULON = SLON + RDIFF 














VION = SLON - RDIFF 












30 


IF (hLON .LT. ULON .AND. 


hLON .GT. 


VLON) 


GO 


TO 


35 



35 


CONTINUE 






C 








C 


DETERMINE 


ASF 


CORRECTOR 


C 










ASF = 0.0 








IF (ID .EQ. 


1) 


ASF = T ABLEW (I, J) 




IF (ID .EQ. 


2) 


ASF = T ABLEX (I, J) 




IF (ID .EQ. 


3) 


ASF = T ABLEY (I, J) 



RETURN 



END 

=== = ==== = == === = = = ===== = = = = === = = = = == = = = = = ^= = = = ^ = ^^=: = =: = = = = 
SUBROUTINE INVERT (A,?INV, GLAT 1 , GLON 1 , GLAT2 , GLON2 ,FAZ, 
1 BkZ,S) 

C 



See Appendix 3 for subroutine INVERT. 



DATA SET EXAMPLE 



PROGRAM LORTAB 



SHIP'S POSITION 
LATITUDE LONGITUDE 



OBSERVED LORAN RATES & ERRORS 



9940-Y 



EY 



9940-W 








(D- 


H-S) 




(D-M 


-s) 


ijxsBc) (;asec) 


(yisec) 


(jisec) 


36 


43 


45. 800 


121 


55 


27. 160 


42789.34 -0.29 


16294.04 


0. 34 


36 


44 


3.400 


121 


55 


32. 340 


42791.13 -0. 18 


16293.46 


0. 30 


36 


44 


21. 180 


121 


55 


37.390 


42793.04 -0.18 


16292. 73 


0. 41 


36 


44 


37.490 


121 


55 


46. 950 


42795.13 -0.38 


16292.03 


0.34 


36 


44 


53. 260 


121 


55 


57. 710 


42796.93 -0.31 


16291.43 


0. 13 



APPENDIX E 



PROG BAM ASFSEL 



C PROGRAM ASFSEL 

C 

C PROGRAM DETERMINES OBSERVED ASF CORRECTORS BY SCANNING 

C DATA AT 1 DEGREE LATITUDE AND LONGITUDE INTERVALS. THE 

C ASF CORRECTORS ARE DETERMINED BY SUBTRACTING THE 

C CALCULATED TD USING THE SEAWATER SECONDARY FACTOR FROM 

C THE OBSERVED TD RATES. 

C 

IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION AXIS ( 13) ,RF (13) 

DIMENSION ASFCR 1 (16,2 6) ,ASFCR2 (16,2 6) , 

1 IN01 (16,26) ,IN02 (16, 26) 

DATA RHOSEC,PI, UNCOV, RN/2 . 06 264 80624 71 DOS , 

1 3. 14 15 9265 35 89 8D0, 29 9. 79 245800,1 .00033800/ 

DATA XMEAN1,XMEAN2,VAR1,VAR2/0. OODO,O.OODO, 

1 0.00D0,00.0D0/ 

DATA XMEAN3,XME AN4/0. OODO,O.OODO/ 

C 

CCCCCC* ***♦***♦♦♦**=>'* ♦♦**♦♦♦♦**#♦* 

DATA AXIS/6. 3782064D0 6,6. 378383 D06 , 6 . 37 7397 1 55D06 , 

1 6.37816 D06, 6. 378 16006,6. 378249145006,6. 378165D06, 

2 6.373166006,6.378165 D06 , 6 . 373 1 4 5D06 , 6 . 377 563 4 D06 , 

3 6.378245006,6.3781350006/ 

DATA RF/6. 3565838 D06, 2.97 D 02 , 2 . 99 152 3 1 2 85D02 , 

1 2.9825002,2. 98 247 16 7 427D 02 , 2 . 9 3 465D02 , 2 . 9 8 25D0 2 , 

2 2.983 D02, 2. 983002,2.9825 D02, 6. 3562569D06 , 

3 2.983 D02, 2.9826D02/ 

C*****ELII?SOID OPTION NUMBER 

C 1. CLARKE 1866 8. MERCURY 
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c 


2. INTERNATIONAL (HAYFORD) 


9. 


MARSHALL ISLAND 


c 


3. BESSEL 1841 


10. 


NAVY 8D 


c 


4. AND (AUSTRALIAN) 


11 . 


AIRY 


c 


5. 1967 REFERENCE 


12. 


KRASSOWSKI 1940 


c 


6. CLARKE 1880 MOD 


13. 


WGS 1972 


c 


7. SAO 







c 

ccc ccc* ♦♦♦♦♦♦ #♦♦♦ * 

C CC1-2 = ELLIPSOID NUMBER (K) 

C CC3-5 = NUMBER OF POINTS ALONG COAST (IREC) 

C 

READ(5,100) K,IREC,ID 1,ID2 
100 FORMAT { IX, 12, 14 ,213) 

TWOPI=2.*PI 
A = AXIS (K) 

F=1./RF (K) 

IF (F. LT.3. D-3) F=(A-1./?)/A 

C 

C READ MASTER AND SLAVE STATIONS POSITIONS 

C THE FIRST RECORD IS THE NUMBER OF MASTER AND SLAVE 

C STATIONS FOR THE PARTICULAR CHAIN. 

C 

WRITS (6,202) 

202 FORMAT(IHI) 

WRITE (6,201) 

201 FORMAT (1H ) 

READ (5, 105) ILATM,IMI NM ,RSECM , I LONM , ILMINM , HLSECH 
WRITE (6,105) ILATM,IMINM,RSECM, ILONM, ILMINM , RLS ECM 
105 FCRMAT(1X,I3, 1X,I2,1X ,?5. 2, 1X,I4, 1X,I2, 1X,?5,2) 

RHAST = { (lABS (ILATM) ♦ 60 + IMINM) * 60 + RSECM) / 

1 RHOSEC 

IF (ILATM .LT. 0) RMAS T = -RMAST 

RMASTL = ( (lABS (ILONM) * 60 + ILMINM) *60 +RLSECM) / 
1 RHOSEC 



IF(ILONfl .GT. 0) RMASTL = TWOPI - RMASTL 
C 

C READ THE FIRST SLAVE STATION AND CHANGE THE LATITUDE 

C AND LONGITUDE INTO RADIANS. 

C 

READ (5, 101) ILATS, IMI NS,RSECS,ILONS,ILMINS, RLSECS , 

1 DELAY1 

WRITE (6,101) ILATS,IMINS,RSECS, ILONS , ILMINS , RLSSCS , 

1 DELAY1 

101 F0RMAT(1X,I3, IX ,12 , 1 X , F5. 2, 1 X, I 4 , 1 X , 1 2, IX, F5 , 2 , F9 , 2) 

RSLAV1 = ( (TABS (ILATS) ♦ 60 + IMINS) * 60 + RSECS) / 

1 RHOSEC 

IF (ILATS .LT. 0) RSLAVl = -RSLA71 

RSLAL1 = ( (TABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 

IF (ILONS .GT. 0) RSLAL1 = TWOPI - RSLALl 
C 

C READ THE SECOND SLAVE STATION AND CHANGE THE LATITUDE 

C AND LONGITUDE INTO RADIANS. 

C 

READ (5, 101) ILATS, IMINS, RSECS, I LONS, ILMINS, RLSECS, 

1 DELAY2 

WRITE (6,101) ILATS, I MINS, RSECS, ILONS, ILMINS, RLSECS, 

1 DELAY2 

RSLAV2 = ( (lABS (ILATS ) * 60 + IMINS) * 60 + RSECS) / 

1 RHOSEC 

IF (ILATS .LT, 0) RSLA V2 = -RSLAV2 

RSLAL2 = ( (TABS (ILONS) * 60 + ILMINS) * 60 + RLSECS) / 
1 RHOSEC 

IF (ILONS .GT. 0) RSLAL2 = TWOPI - RSLAL2 

WRITE (6,201) 

C 

C READ THE RECORDS POSITIONS AND OBSERVED RATES, THEN 

C COMPUTE THE TOTAL DISTANCE. 
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c 

DO 550 I=1,IREC,1 

READ (5, 104) ILAT,IMIN , RSEC , ILON , ILMI N , RLSEC , 

1 RATE1,RATE2 

104 FORMAT (IX, 13,13 ,F6. 3, 15 ,13 , F6 . 3 , 3 X, 2F1 0. 2) 

C 

C CONVERT LAT AND LONG TO RADIANS. 

C 

RPOST = ((IA3S(ILAT) ♦ 60 + IMIN) * 60 + RSEC) / 

1 RHOSEC 

IF(ILAT .LT. 0) RPOST = -RPOST 

RPOSTL = ((lABS(ILON) ♦ 60 + ILMIN) ♦ 60 + RLSSC) / 

1 RHOSEC 

I? (ILON .GT. 0) RPOSTL = TWOPI - RPOSTL 
C 

C COMPOTE DISTANCES AND AZIMUTHS FROM THE OBSERVED POINT 

C 

C MASTER 

C 

CALL IN7ER1 (A, F, RPOST , RPO STL , RM AS T , R M ASTL, F AZM , B AZM , 

1 DISTM) 

ODISTM = (RN * DISTM) / DNCOV 
CALL SSCF AC (ODISTM, SFM) 

C 

C SLA7E1 

C 

CALL INVERT ( A , F , RPOST , RPO STL, RS L A V 1 , RSL AL 1 , F AZ 1 ,3AZ1, 

1 DIST1) 

UDIST1 = (RN * DIST1) / ONCOV 
CALL SECFAC (UDI ST1 ,SF 1) 

C 

C SLAVE2 

C 

CALL INVERT (A,F , RPOST , RPO STL, RS L A 72 , RSLAL2, F AZ2 , 3 AZ 2, 
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1 DIST2) 

UDIST2 = (RN * DIST2) / UNC07 
CALL SECFAC (UDIST2,SF 2) 

C 

C DETERMINE WHICH LAT AND LONG THE ASF CORRECTOR 

C IS ASSIGNED TO. 

C 

CALL ASSIGN (RPOST,RPO STL, JN1, JN2) 

C 

C COMPUTE THE RATES AND COMPARE TO THE OBSERVED RATES 

C 

TDM1 = 0DIST1 - ODIST M + SF1 - SFM + DELAYl 
DIFF1 = TDM1 - RATE1 
C 

TDM2 = ODIST2 - ODIST M + SF2 - SFK + DELAY2 
DIFF2 = TDM2 - RATE2 
C 

C SUM THE DIFFERENCES TO THE MATRIX AND COUNT THE NUMBER 

C OF ASF CORRECTORS FOR EACH BLOCK TO LATTER DETERMINE 

C THE MEAN. 

C 

ASFCR1 (JN1, JN2) = ASF CR 1 ( JN 1 , JN 2) + DIFF1 
ASFCR2( JN1, JN2) = ASF CR2 ( JN 1 , JN 2) + DIFF2 
IN01 (JN 1, JN2) = IN01 ( JN1, JN2) + 1 
550 IN02 (JN1, JN2) = IN02 ( JN1, JN2) + 1 
C 

C DETERMINE THE MEAN ASF CORRECTOR FOR EACH LAT AND LONG 

C 

DO 650 1=1 , 16,1 
WRITE (6,201) 

DO 600 J=1,26,1 

IF(IN01(I,J) .EQ. 0) GO TO 580 
ASFCR1(I,J) = ASFCai(I,J) / IN01(I,J) 

IF(IN02(I,J) .EQ. 0) GO TO 590 
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580 



ASFCR2(I,J) = ASFCR2(I,J) / IN02(I,J) 

590 IF (ASFCR1 (I, J) .EQ. 0.0) ASPCR1(I,J) = 9.99 

600 IF (ASFCR2 (I, J) .EQ. 0.0) ASFCR2(I,J) = 9.99 

650 CONTINUE 
C 

C WRITE THE CORRECTORS IN MATRIX FORMAT 

C 

c 

IF(ID1 .EQ. 1) WRITE (6, 1000) 

IF(ID1 .EQ. 2) WRITE (6, 100 1) 

IF(ID1 .EQ. 3) WRITE ( 6, 1002) 

WRITE (6,201) 

C 

WRITE (6,300) ( (ASFCR1 (I,J) ,1 = 1,16) , J=1,26) 

DO 700 L=1 , 5, 1 
WRITE (6,201) 

201 FORHAT(1H ) 

700 CONTINUE 
C 

IF(ID2 .EQ. 1) WRITE ( 6, 1000) 

IF(ID2 .EQ. 2) WRITE (6 , 100 1) 

IF(ID2 .EQ. 3) WRITE ( 6, 1002) 

WRITE (6,201) 

WRITE(6,300) ( (ASFCR2 (I,J) ,1 = 1,16) ,J = 1,26) 

300 FORMAT(16F6.2) 

1000 FORMATC TABLE FOR 99 40-W ’) 

100 1 FORMATC TABLE FOR 99 40-X ') 

1002 FORMATC TABLE FOR 99 40-1 ») 

STOP 

END 

======= = ======= = == ===== = ========== = ==== = ====== 

SUBROUTINE SECF AC ( UTD 1ST, 3 F) 

C 

C 
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C THIS ROOTINE WILL COMPUTE THE SEA SECONDARY FACTOR 

C UTDIST = TOTAL DISTANCE 

C SF = SECONDARY FACTOR 

C 

IMPLICIT REALMS (A-H, 0-Z) 

C 

C 

C COEFFICIENTS 

C 



AO 


= 


129.04323 


A1 




-0.40758 


A2 


= 


0.00064576813 


BO 


= 


2.741282 


B1 


= 


-0.011402 


B2 


= 


.000327748 15 



C 

IF (UTDIST .GT. 537) GO TO 10 

SF = ( BO / UTDIST) + B 1 + ( B2 * UTDIST) 

GO TO 20 

10 SF = ( AO / UTDIST) ♦ A1 + ( A2 * UTDIST) 

C 

20 CONTINUE 

RETURN 
END 

C= ============ === = == ========== ====== = = = ===== === = = === = 

SUBROUTINE ASSI GN (RL A T, RLON, I , J) 

C 

C 

C SUBROUTINE ASSIGN SELECTS THE COLUMN AND ROW FOR 

C THE LATITUDE AND LONGITUDE OF THE RECORD. 

C RLAT = POSITION LATITUDE IN SECONDS 

C RLON = POSITION LONGITUDE IN SECONDS 

C I = COLUMN 



C J = ROW 

C 

IMPLICIT REAL+8 (A-H^O-Z) 

DATA RHOSEC,PI/2.062648062471D05, 3. 1 4 1 59265 35898D0/ 
C 

TWOPI = PI * 2. 0 

C CONVER RLAT AND RLON TO SECONDS 

C 

HLAT = RLAT * RHOSEC 
HLON = TWOPI - RLON 
HLON = HLON * RHOSEC 
C 

C STARTING LAT AND LONG FOR SEARCH 

C LAT = 37/05/00.0 LONG = 122/05/00.0 

C 

SLAT = 133200.0 
SLON = 439500.0 
C 

C LATITUDE 

C 

J = 0 

DO 10 IC = 1,26 ,1 
SLAT = SLAT - 60.0 
RULAT = SLAT + 30.0 

RLLAT = SLAT - 30.0 

10 IF (HLAT .LT. RULAT .AND. HLAT .GE. RLLAT) GO TO 15 

15 J = IC 

C 

C LONGITUDE 

C 

1=0 

DO 20 IC = 1,16,1 
SLON = SLON - 6 0.0 
RLLON = SLON + 30.0 
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RRLON = SLON - 30.0 

20 IF(HLON ,LT. RLLON .AND. HLON .GS. RRLON) GO TO 25 
25 I = IC 
RETURN 
END 

C ================== ======== ========== = =========== ========== 

SUBROUTINE INVERT ( A, F INV, GLAT 1 , GLON 1 , GLAT2 , GLON2 , FAZ , 
1 BAZ,S) 

C 



See Appendix B for subrcurine INVERT. 
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DATA SET EXAMPLE 



PROGRAM ASFSEL 



TABLE FOR 9940-Y 
50* 



12 20 / 00 ' 



-0.53 

45* -0.56 -0.51 -0.46 -0.40 

-0.39 -0. 54 -0.39 -0.55 - 
-0. 63 -0.57 - 
-0.54 - 



40* 



35' 

360 



0.41 

0.44 

0.45 

0.49 

0.50 

0.62 
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TABLE FOR 9940-Y 

1210 / 551 1210 / 50* 

50 * 





1 

0 

• 

1 

0 

• 

1 

0 

• 


-0.37 -0.60 




45* 


-0.51 -0.50 


-0. 4 3 


-0. 44 




-0.38 -0.49 -0.50 


-0.46 -0.52 


-0. 46 




-0.54 -0.51 -0.49 








-0.54 -0.57 -0.49 








-0.55 -0.55 -0.50 






40* 


-0.67 -0.61 -0.50 








-0.64 -0.65 







35* 



0.44 
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